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Abstract 

In this thesis we consider the phenomenology of supergravity, and in particular the 
particle called "gravitino". We begin with an introductory part, where we discuss 
the theories of inflation, supersymmetry and supergravity Gravitino production is 
then investigated into details, by considering the research papers here included. 

First we study the scattering of massive W bosons in the thermal bath of par- 
ticles, during the period of reheating. We show that the process generates in the 
cross section non trivial contributions, which eventually lead to unitarity breaking 
above a certain scale. This happens because, in the annihilation diagram, the lon- 
gitudinal degrees of freedom in the propagator of the gauge bosons disappear from 
the amplitude, by virtue of the supergravity vertex. Accordingly, the longitudinal 
polarizations of the on-shell W become strongly interacting in the high energy limit. 
By studying the process with both gauge and mass eigenstates, it is shown that 
the inclusion of diagrams with off-shell scalars of the MSSM does not cancel the 
divergences. 

Next, we approach cosmology more closely, and study the decay of a scalar 
fleld S into gravitinos at the end of inflation. Once its mass is comparable to 
the Hubble rate, the fleld starts coherent oscillations about the minimum of its 
potential and decays pertubatively. We embed in a model of gauge mediation 
with metastable vacua, where the hidden sector is of the O'Raifeartaigh type. First 
we discuss the dynamics of the fleld in the expanding background, then radiative 
corrections to the scalar potential V{S) and to the Kahler potential are calculated. 
Constraints on the reheating temperature are accordingly obtained, by demanding 
that the gravitinos thus produced provide with the observed Dark Matter density. 
We modify consistently former results in the literature, and find that the gravitino 
number density and Tr are extremely sensitive to the parameters of the model. This 
means that it is easy to account for gravitino Dark Matter with an arbitrarily low 
reheating temperature. 
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Chapter 1 

Inflationary cosmology 



"Tyger! Tyger! Burning bright 
In the forests of the night, 
What immortal hand or eye 
Could frame thy fearful symmetry?" 

- William Blake, "The Tiger" 

"E quindi uscimmo a riveder le stelle. " 

- Dante Alighieri; Inferno XXXIV, 139 

Cosmology, (from the Greek KoafxoXojia, namely -Koaiioq, kosmos, "Uni- 
verse" and -XojLa, -logia, "study") has extremely ancient origins, which can be 

dated back to the prehistoric era. Though the word Cosmology is recent (first used 
in 1730 in Christian Wolff's Cosmologia Generalis), the study of the Universe has a 
long and complex history. It involves science, religion, philosophy and esoterism. 

Thousands years ago the Babylonians were skilled in Astronomy, building up a 
tradition that was further developed by the Ancient Greeks. The latter were the first 
to build a cosmological model within which to interpret the motions of planets and 
stars. In the fourth century BC, they believed that the stars were fixed on a celestial 
sphere which rotated about the spherical Earth every 24 hours. The planets, the 
Sun and the Moon, would have then moved in the ether between the Earth and the 
stars. 

This model was refined, and in the second century AD it lead to Ptolemy's sys- 
tem. It was based on religious and philosophical fundaments, as the circular motion 
of stars and planets was justified by the belief that perfect motion should be in cir- 
cles. The earth was put at the centre of the Universe, by choosing a reference frame 
which turned out not to be practical. Four geometric devices, the deferent, epicy- 
cle, eccentric and equant had to be introduced to describe the planetary motions 
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correctly. The resulting system was extremely complicated. 

However, Ptolemy's model was successful in describing the celestial motions 
with great accuracy. Much later, in the 16th century, Copernicus assumed that the 
Earth was not at the centre of the Universe, but rotated and moved in a circular 
orbit about the Sun. Unfortunately, his model could not match the accuracy of 
Ptholemy's system, which was still regarded as the standard cosmological model. 

Nevertheless, after a few decades Galileo discovered that there are moons or- 
biting the planet Jupiter. This was a first observational evidence against the Earth- 
centred system. The final strike of genius came from Johannes Kepler. Kepler 
postulated that the planets moved in ellipses, not perfect circles, about the Sun. 
The Heliocentric model soon replaced Ptolemy's system, but it still was an entirely 
phenomenological model. The theoretical justification for Kepler's laws came after 
about one century, when Newton showed that elliptical motion could be explained 
by his law for the gravitational force. 

After these fundamental discoveries, the interest in the solar system was then 
pushed forward to larger scales. The philosopher Immanuel Kant, among the others, 
proposed that the Milky Way was just one "island Universe" or galaxy, and that 
beyond it must be other galaxies. 

Since the formulation of the theory of General Relativity by Einstein in 1915, 
it became possible to discuss the evolution of the Universe from a fully consis- 
tent physical viewpoint. The Russian mathematician and meteorologist Alexander 
Friedman found in 1917 that the Einstein equations could describe an expanding 
Universe. This implied an instantaneous birth of the Universe, about ten thousand 
million years ago in the past, with creation of all the matter at just one instant. 
The British astronomer Fred Hoyle called it the " Big Bang" . 

In 1929, Edwin Hubble discovered through observations of the redshift of galax- 
ies [3] that the Universe is expanding. George Gamow and his collaborators found 
in 1946 that nucleosynthesis requires a hot and extremely dense state at the origin 
of the Universe. This statement is now regarded as the standard Big Bang scenario. 
Remarkably, they also predicted that the actual Universe is filled with a background 
radiation of the black body-type. This was confirmed by the experimental discov- 
ery of the Cosmic Microwave Background Radiation (CMBR or simply CMB) by 
Penzias and Wilson in 1965 [4j. 

In the standard Big Bang scenario, the state of the Universe (either radiation- 
or matter-dominated) consists of a decelerated expansion, due to the negative sec- 
ond derivative of the scale factor. However, several cosmological problems such as 
flatness and horizon problem plague this description. The theory of inflation was 
originally proposed in 1980 by Alan Guth in [5] and by Katsuhiko Sato in [6J, as a 
mechanism for resolving these problems. Contemporary with Guth, Alexei Starobin- 
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sky argued that quantum corrections to gravity would replace the initial singularity 
of the Universe with an exponentially expanding state [7j. Einhorn and Sato pub- 
lished a model similar to Guth's and showed that it would resolve the puzzle of 
the magnetic monopole abundance in Grand Unified Theories |8j. Like Guth, they 
concluded that such a model (which is now called old inflation) not only required 
fine tuning of the cosmological constant, but also would very likely lead to a much 
too granular Universe, namely to large density variations resulting from bubble wall 
collisions. 

The bubble collision problem was solved in 1982 by Andrei Linde in [9] and 
independently by Andreas Albrecht and Paul Steinhardt in [10] , in a revised version 
that is named new inflation. 

The basic idea is slow-roll inflation, where instead of tunneling out of a false 
vacuum state, inflation occurred by a scalar field rolling down a potential energy 
hill. When the field rolls very slowly compared to the expansion of the Universe, 
inflation occurs and eventually ends when the potential becomes steeper. However, 
there was still a problem of fine tuning also in this model, since the scalar field would 
need to spend enough time in the false vacuum to provide with a sufficient amount 
of inflation. Eventually, in 1983 Linde [Hj introduced a variant version of the slow- 
roll in which initial conditions of the scalar fields are chaotic. This corresponds 
to chaotic inflation. It predicts that our homogeneous and isotropic Universe may 
be produced in the regions where inflation occurs sufficiently. The basic difference 
between chaotic and old (or new) inflation, is that the latter requires a Universe 
in a state of thermal equilibrium from the beginning, whereas the former does not 
necessarily need this assumption. Moreover, chaotic inflation powerfully solves also 
the problem of initial conditions, since it can start at Planckian densities. 

Since the definition of cosmic inflation, many different models have been con- 
structed [12]. In particular, the trend is now to incorporate particle physics in the 
model, in order to create a fully consistent physical theory [13l [14]. In fact, the 
detailed particle physics mechanism responsible for inflation still needs to be under- 
stood. The issues of preheating, reheating, the generation of density perturbations 
and the creation of structures, involve a whole class of hard-core physics, from su- 
pergravity to string theory. Also the origin and the nature of the inflaton field are 
still unclear and problematic, as in chaotic inflation it is just a scalar particle with 
no defined structure. Attempts to embed this field within a particle physics theory 
are object of intense discussion nowadays. For instance, there are models where the 
inflaton is defined by linear combinations of flat directions in the MSSM [T5] . 

The inflationary paradigm not only solves the flatness and horizon problems, 
but it also generates density perturbations which then concur to the formation of 
cosmological structures in the Universe. This is very important for observations, as 
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it was first noticed in the early 80s by Stephen Hawking [16], Starobinsky [T7] and 
several others [18]. Nearly scale-invariant spectra of cosmo logical perturbations are 
generated by quantum fiuctuations of the infiaton, whose scales are frozen during 
the expansion. Long after the ending of infiation, the scales re-enter the Hubble 
radius. Accordingly, the perturbations which originated during infiation can create 
the structures of the Universe which correspond to a large scale. Such scale-invariant 
spectra were observed by the COBE satellite in 1992 [19], and the recent WMAP5 
data [20] confirm this interpretation. 

In this thesis we cover some aspects of this fascinating, yet problematic subject. 
We shall attempt to describe a number of cosmological issues from the point of view 
of particle physics. In the introductory part we provide with some basics which are 
relevant to our discussion, which then focuses on the research papers here included. 

In Chapter 1 we discuss the theory of cosmological infiation, that provides 
the background for the physical mechanisms which are studied in this work. After 
considering the standard Big Bang cosmology, we review old and new infiation, 
chaotic infiation and reheating. 

Chapter 2 introduces supersymmetry (SUSY) and supergravity (SUGRA). We 
discuss the spontaneous breaking of global and local SUSY, and the gravitino as the 
gauge field of supergravity. We obtain from the general SUGRA Lagrangian the 
interactions of the gravitino. 

Chapter 3 refers to the paper [Ij. We investigate the scattering of W bosons 
into gravitino and gaugino in the broken phase, by using both gauge and mass 
eigenstates. Differently from what is obtained for unbroken gauge symmetry, we 
find in the scattering amplitudes new structures, which can lead to violation of 
unitarity above a certain scale. We show that the longitudinal polarizations of the 
on-shell W become strongly interacting at high energies, and show that the inclusion 
of diagrams with off-shell scalars of the MSSM does not cancel the divergences. 

Chapter 4 considers the production of Dark Matter gravitinos via the decay of a 
scalar field at the end of infiation. We discuss the dynamics of the field and the model 
of gauge mediation in which it is embedded. It is shown that the gravitino density 
is extremely sensitive to the parameters of the hidden sector. For the case of an 
O'Raifeartaigh hidden sector, the observed Dark Matter density can be explained by 
gravitinos even for low reheating temperatures Tr 10 GeV. This chapter examines 
the contents and results of the paper [2]. 

In the Appendix we summarise our conventions for spinors, propagators and 
spin sums. We provide also with a list of the Feynman rules for the MSSM and for 
supergravity which are relevant for this thesis. 



1.1. THE HISTORY OF THE UNIVERSE IN A NUTSHELL 5 

1.1 The history of the Universe in a nutshell 




Figure 1.1: Schematic history of the Universe. 



Since the beginning of time (provided it did have a beginning), the Universe 
has passed through several epochs, which are summarised in Fig 11.11 and hsted 
chronologically below. Some of them are described by established physics that is 
verified by experimental data, others by hypotheses and conjectures based on physics 
beyond the Standard Model. Some are still completely obscure to our understanding 

m- 

• Nobody actually knows anything certain about the beginning. Was it a period 
that can be described by quantum gravity? What about time, did it exist 
or not? These are still open questions, which belong to the so-called pre- 
inflationary cosmology. 

• Inflation, namely a period of exponential expansion of the Universe. This 
model is nowadays widely accepted by the scientific community, since it cures 
the problems of the standard big bang cosmology in agreement with the ob- 
servational data. The quantum fluctuations during inflation were indeed im- 
printed on the metric. They can be observed as CBR fluctuations in the power 
spectrum, namely as deviations from homogeneity and isotropy in the matter 
distribution. 



CHAPTER 1. INFLATIONARY COSMOLOGY 



End of inflation, particle production. Once that the exponential expansion has 
come to an end, the inflaton field releases the energy of the Universe. Photons 
and other elementary particles are produced via several mechanisms, which 
are addressed in the next sections. 

Baryogenesis. The confinement of quarks takes place and baryons are formed. 
The excess of matter with respect to antimatter is created. 

The Universe reaches thermal equihbrium and it is cooled down adiabatically. 

Several phase transitions take place in this epoch: breaking of grand unified 
theories (GUT), electroweak (EW) symmetry, supersymmetry (SUSY), phase 
transition from free quark-gluon phase to confinement in quantum chromody- 
namics (QCD) and so on. Possible formation of topological defects. 

Decoupling of neutrinos from the cosmological plasma, when the Universe is 
1 sec old and T ^ 1 MeV. 

Big Bang Nucleosynthesis (BBN) is the epoch when the light elements ^H, 
^He, "^He, Li were formed. It occurs in the time interval from 1 s to 200 s, 
and T 1 — 0.07 MeV. It is probably the only mechanism in the early Universe 
for which we have a very good agreement of theory with data. 

Structure formation starts when the cosmological matter turns from relativistic 
to non-relativistic. This corresponds to the epoch of equivalence between 
radiation and matter domination. It takes place at T ~ eV, with redshift 

Decoupling from matter of the Cosmic Microwave Background (CMB), which 
after that has propagated almost freely in the Universe. This is called the 
period of hydrogen recombination, aiT 0.2 eV or Zrec ~ 10^. Baryons began 
to fall into already evolved seeds of structures generated by Dark Matter. 

A variegated cosmological jungle of stars, planets, black holes, worm holes 
and other cosmological objects is formed. The present time corresponds to 
t ^ 12 - 15 X lO^yr and T = 2.3 x lO-^^Qgy = 2.7 K. The redshift z is 
defined by 2; + 1 = a{to)/a{t), where a(to) is the scale factor of the universe at 
the actual time to- Accordingly, zq = z{to) = 0. 
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1.2 The standard Big Bang cosmology 



Let us assume the cosmological principle [T2]: at large scales, the Universe is homo- 
geneous and isotropic. Technically speaking, this is achieved by using the Friedman- 
Robert son- Walker (FRW) metric 



ds'^ = g^udx'^dx'^ = —di^ + a^{t) 



dr^ 



1 — fcr^ 



+ r\d9^ + sm^ edcj)^) 



where a{t) is the scale factor, which measures the radius of the Universe according 
to the flow of the cosmic time t. The constant k is the cosmological curvature, which 
can take the values (fiat Universe), 1 (closed of spherical Universe) or -1 (open of 
hyperbolic Universe). 

Besides the geometrical structure of the Universe, the matter content is pri- 
marily important for its evolution. This is defined by the equation of state between 
the energy density p{t) and the pressure p(t). Namely, for a radiation-dominated 
Universe we have p = p/3, whereas if it is matter (or dust) dominated, p = 0. Given 
the metric and the matter content, the dynamics of the Universe is determined 
by solutions of the Einstein equations in General Relativity. With a cosmological 
constant^, they read as 



SttGT, 



'1.21 



Here R^y is the Ricci tensor, R is the Ricci (curvature) scalar, T^^ is the energy 
momentum tensor and Gfj_y is the Einstein tensor. G is Newton's gravitational 
constant and it is related to the Planck mass mp = 1.2211 x 10^^ GeV through 
mp = {hc/Gy/'^ (or either mp = (l/G)^/^ in natural units). Throughout this thesis, 
we will use extensively also the reduced Planck mass, defined as Mp = mp/\/87i = 
2.43 X lO^s GeV. 

Remarkably, the Einstein equations (11. 2p summarise the main challenges of 
contemporary cosmology. The left hand side contains the geometric structure of the 
Universe, which can be modified by either using the cosmological constant or by 
properly modifying the Einstein tensor. The latter approach is discussed in theories 
of modified gravity [22], which have been recently proposed as an alternative to dark 
energy. On the other hand, the energy momentum tensor T^y encodes the matter 
constituents, among which Dark Matter. Accordingly, the topics of this thesis deal 
with the right hand side of (11. 2p . 



"'^The cosmological constant A was originally postulated by Einstein to achieve a stationary 
Universe. Even though it is now clear that our Universe is expanding (Einstein called it the 
"biggest blunder" of his life), since A is directly related to dark energy, this was clearly another 
strike of genius of his. 
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By using the FRW metric (11. ip . and ignoring the cosmological constant, the 
Einstein equations become the two Friedman equations. From the 00 component of 
the Eistein tensor we obtain 



u2 ^ n Q^ 



where the Hubble expansion rate is defined as H = a/ a. The components G^^ 

Q22 ^ ^33 

give the dynamics of the scale factor a{t): 



4:71 

— ip + Sp). (1.4) 



a Snip 

In the above, we assume the perfect fiuid form for the energy- momentum tensor 

T^"" = {p + p)u^'u'' + pg^"" . (1.5) 

By combining (11. 3p and (II. 4p . one finds the continuity equation, 

p + 3R{p + p) = 0, (1.6) 

which is the General Relativity version of energy-momentum conservation. The 
Friedman equation (11.31) can also be rewritten as 

^^-1 = ^, where n^-^ = -^p. (1.7) 

The dimensionless parameter Q measures the amount of deviation of the actual 
energy density p of the Universe from the critical density pc = SH^mp/Svr, that 
corresponds to a spatially fiat geometry. In fact, if = then p = Pc and f2 = 1. In 
this case, the solutions of the Friedman equations are 

a oc p (X a , (1-8) 

for radiation domination and 

a oc t^/^ , p oc (1.9) 



for matter domination. In both cases, this corresponds to an expanding Universe 
with a decelerating expansion, since the second time derivative of the scale factor is 
negative (a < 0) by virtue of Eq. (ll.4p . 
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1.2.1 Problems of the standard scenario 

The cosmological models (also named "Friedman models") above defined are consis- 
tent from a theoretical viewpoint, but they are plagued by a number of phenomeno- 
logical problems [T3] . 

Let us discuss them in this section. The flatness problem originates from 
Eg. (11. 71) . Since the term a^H^ always decreases, the matter density Q shifts away 
from 1 with the expansion of the Universe. In contrast, observations seem to suggest 
that nowadays Q ~ 120] ■ Accordingly, it must have been very close to 1 in 

the past. It can be shown that in order to have 0.1 ^ fl ^ 2, consistently with 

2 

the experimental data, the early Universe must have had Ifi — 11 ^ 10~^^7^. This 
implies that at the Planck epoch when T mp. 



\n 



^-1 

pc 



10-^^ (1.10) 



Thus if the initial density would have been larger than pc by e.g. a factor of 10"^^, 
the Universe would have collapsed long ago. Conversely, if the density were smaller 
than Pc by the same factor, the expansion would have been so rapid to forbid the 
formation of structures. Therefore we are left with an extreme fine-tuning of the 
initial conditions, namely with a fairly unnatural coincidence in the history of the 
Universe. 

The horizon problem deals with causality. First consider a comoving wave- 
length A and a physical wavelength aX such that aX ^ (namely, it is inside the 
Hubble radius H^^). Since in the Friedman models a (x with < p < 1, we 
find that aX oc t^. However, H~^ oc t and the physical wavelength becomes much 
smaller than the Hubble radius, defining a causal region which is very small within 
the horizon. 

More into detail, with the notations of Ref. j23j, the particle horizon Dii{t) is 
defined as the region where the light travels from the moment of the big bang , 

/■* dr 

Duit) = aifyduit) , dH{t)= (1.11) 

Jt, «(^) 

where dnit) is the comoving distance. The photons we observe in the CMB were 
emitted at the time of decoupling, and D nit dec) = o-ii dec) dn (idee) corresponds to the 
region where photons have interacted at that time. If the comoving distance at the 
actual time (to) is dnito), we have 

dnitdec) ftoY^^^ fm^Y^^ , 

HKdec) oc — -r ^ 2 X 10"^ . 1.12 

duito) \tdec VlQlV 
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This implies that the regions where there is causal interaction of photons are small. 
However, in the CMB we actually see photons which thermalise to the same temper- 
ature in all regions. This is in contradiction with the standard big bang cosmology. 

The origin of large scale structure is also problematic. The anisotropics in 
the last scattering surface (LSS) have been observed by the COBE satellite. Their 
amplitudes are small and almost scale-invariant. In the standard cosmology, it is 
impossible to generate them between the big bang and the LSS, because they spread 
too fast and spoil the formation of structures. 

Finally, the monopole problem regards the creation of many unwanted relics 
(such as monopoles, topological defects and cosmic strings) due to supersymmetry 
breaking [13]. The cosmic string, in particular, predicts the primordial formation 
of gravitinos and moduli fields, whose energy density decreases like matter, namely 
as a~^. The radiation energy density decreases as in the radiation-dominated 
era, therefore such relics would be dominant in the Universe. This is of course 
inconsistent with observations. 

1.3 The inflationary model 

The physical mechanisms which are analysed in Chapters |3] and |4] take place right 
after the end of inflation, which therefore provides a common background. Following 
Refs.|13[ |23] . we here provide with an introductory description of the inflationary 
paradigm. 

The basic idea is simple, yet extremely powerful. Since the problems of the 
standard Big Bang cosmology are related to an always decelerating expansion, let 
us assume that in the early Universe there has been an early stage with accelerated 
expansion. In other words, we impose 



Eg. (11. 131) means that a (and therefore aH) increases during inflation. Accordingly, 
the comoving Hubble radius (aH)~^ now decreases, and the a^H^ term in 



a > , 




that by virtue of (ll.4p implies 



p+3p<0. 




n - 1 



k 




a 



2H2 ' 



increases during inflation and Q becomes rapidly 0{1). After inflation ends, Q starts 
to decrease like in the standard scenario. However, if inflation lasts long enough, Q 
remains close to unity until the present day. The flatness problem is now avoided. 
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Moreover, during inflation the scale factor evolves as a oc t^, with p > 1. This 
implies that the physical wavelength a\ is pushed outside the Hubble radius, which 
grows linearly with time. Namely, causality is now guaranteed in regions which 
are much larger than the horizon, and the horizon problem is solved. Also in this 
case, however, inflation should last for enough time, since after inflation the Hubble 
radius starts to grow again faster than the physical wavelength. Since we need that 
before decoupling the photons cover a comoving distance that is much larger than 
the distance after the decoupling, the following condition must be satisfied. 



This happens if the Universe expands times during inflation [T2]. 

Large scale structures can be generated after inflation, since the comoving 
Hubble radius decreases during the inflationary expansion. This means that the 
nearly scale-invariant perturbations, which are needed for the creation of cosmo- 
logical structures, are causally related and small quantum fluctuations are thus 
generated. After the scale is pushed outside the Hubble radius during inflation, the 
perturbations can be described as classical. When inflation ends, the evolution of 
the Universe is described by the standard Big Bang model, and the comoving Hubble 
radius begins to increase. Then the scales of perturbations cross inside the horizon, 
and causality follows. The small perturbations originated during inflation appear 
as large-scale perturbations after this second horizon crossing. This produces the 
seeds of density perturbations which are observed in the anisotropics of the CMB. 

We conclude with a comment on particle production. During the inflationary 
epoch the energy density decreases very slowly. Indeed, a oc t^, where p > 1, 
implies H oc oc a~^^^, and p oc a~^/^. At the same time, the energy density of 
massive particles decreases much faster (oc a~^), thus the particles are red-shifted 
away during inflation and diluted. Also the monopole problem is thus solved. 



1.3.1 Slow-roll conditions and inflationary models 

After considering the main idea of inflation and its cosmological effects, we now 
discuss the theory more closely. Let us consider a homogeneous scalar field if, 
which we will call the inflaton. The potential V{ip) makes the Universe expand 
exponentially. The energy density and the pressure density of this particle are 
written as: 





p = -0' + V{ip) , 



(1.17) 
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Substituting this into fll.Sp and (11. 6p . we get 

'1 



H2 



3mp 



+ Vi^p) 



:i.i8) 



and 

+ 3Hv9 + V'(v^) = , (1.19) 

where we have neglected the curvature term k/a"^. The above means that during 
inflation, Eq. (ll.l4|) provides with ip"^ < V{(f), namely the potential energy of the 
inflaton dominates over its kinetic energy. Accordingly, since we need a very flat 
potential to guarantee a sufficient amount of inflation, we impose the following slow- 
roll conditions: 

vi'«V(¥.), (1.20) 
|<^|<|3H0|. (1.21) 

By introducing the fundamental slow-roll parameters, 

'V'V m%V" 



m. 



1671 \ V ' ' Svr 1/ 



1.22) 



it can be shown that Eqs. (ll.20p and (11.211) imply e ^ 1, |?7| <^ 1. In this limit, 
Eqs. (ll.l8|) and (I1.19P can be rewritten respectively as follows, 

H^^^n^^), (1.23) 
3Hv9 ^ -V\^) . (1.24) 

The above are called the slow-roll equations. Accordingly, inflation ends when either 
e or ry grow enough so that they approach order unity. We remark that the conditions 
e ^ 1 and jryl ^ 1 are just constraints on the shape of the potential, and that they 
do not necessarily imply the slow-roll equations. 

To measure the amount of inflation, the standard quantity is the number of 
e-foldings, defined as 

N = \n'^= j mt, (1.25) 
fli Ju 

where the subscripts i and / denote respectively the quantities at the beginning and 
at the end of inflation. In order to solve the flatness problem, — 1| ^ 10"^*^ right 
after the end of inflation. The ratio of this quantity between the initial and final 
phase of inflation is 



< 



— 1| \af 



e-''' , (1.26) 
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if one assumes that H is nearly constant in that period. From the above equation, if 
— 1| ~ 1, we find that the number of e- foldings has to be ^ 70 to solve both 
the flatness problem and the horizon problem. 

In the next subsections we will sketch some inflationary models, which differ 
by the particle content and the scalar potential V{(p). In particular, we will consider 
chaotic and hybrid inflation. There exists a number of other models, for instance 
natural inflation [23], but we will not describe them here. 



Chaotic inflation 



This model was proposed in 1983 by Andrei Linde ^llj. It is defined by chaotically 
distributed initial conditions and by a potential that can be either quadratic. 



or quartic. 



Vi^) = imV' , (1.27) 

V{^) = \X^\ (1.28) 

thus provided with a self-interaction term. With a quadratic potential, the slow-roll 
equations f ll.23p and f ll.24p can be rewritten as 



H^^ 3H0 + mV~O, (1.29) 

3mp 



which provide the following solutions. 



2V37r 

and 



a tti exp 



mmp 

^^^^-—^t, (1.30) 



:i.3i) 



TT m ( mmp 
tyipi- t 



3mp \ 4A/37r 

(fi is the initial value of the inflaton. What happens physically, is that the inflaton 
finds itself displaced from the true vacuum, to which it rolls back, as in Fig. (1.2). In 
the case of a quartic potential, when > \~^^^mp, Lp has a greater energy density 
than the Planck density, and classical physics cannot describe the evolution of the 
Universe. Instead, if mp/3 < < X~^/'^mp, the values of the field 99 slowly decrease 
and the Universe then inflates. Inflation takes place while the inflaton is displaced. 
The evolution of the scale factor implies the required exponential expansion during 
the initial stages of inflation. Then it slows down since {—mmp / A\/^)t^ grows with 
time. The slow-roll parameters in this case are identical, that is 
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Figure 1.2: Chaotic inflation, limp/3 < 
if < X~^^^mp, the values of the field ip 
slowly decrease and the Universe then in- 
flates. Infiation takes place while the in- 
flaton is displaced, as in figure ^13j . 




Accordingly, inflation ends when ^ rup/An. The initial value of the inflaton field 
is obtained by demanding sufficient inflation, i.e. that ^ 70. This is fulfilled by 
(Pi ^ Smp, since the number of e-foldings is given by 




(1.33) 



The observations of the density perturbations by the COBE satellite [19] provide 
with constraints on the inflaton mass, 

m ^ 10"W , (1-34) 
and on the self-coupling [12], which is very small A ~ 10~^^. 



Hybrid infiation 

Instead of considering only a single inflaton field, one might want to define multiple 
scalar fields [25]. In this case, the potential can be written as 

^ = ^ (x^ - ^) ' + + . (1.35) 

If if"^ is large, this can be approximated by the single-field potential 

^=4Y + 2^V, (1-36) 

since (p tends to roll down toward the potential minimum at x = 0. Inflation does not 
end for the approximated potential fll.36p . however the mass of x becomes negative 
for (fi < = M/g. This implies that the field rolls down to one of the true minima 
at (f = and x = ±M/V^, as illustrated in Fig.(1.3). Inflation ends after the 
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Figure 1.3: Hybrid inflation. The inflaton 
rofls down to one of the true minima at 
if = and X = ±M/\/A. Inflation ends 
after the symmetry breaking for (p < (p^, 
due to the roUing of the field x- 




symmetry breaking for < yjc, due to the rolhng of the field x- ^ the initial value 
of the infiaton is <^i, the number of e-foldings is 



Right after the end of infiation, which cooled the Universe via the expansion, a 
reheating period occurs. The Universe is thermalised, as the potential energy of the 
inflaton is transferred to radiation. In the original idea, the "old reheating" [26] , 
the inflaton decays perturbatively. However, this is not an efficient mechanism for 
baryogenesis at the GUT scale. It was later found that a consistent scenario should 
include a nonperturbative stage called preheating, where an explosive production of 
particles occurs in the early stages of reheating [271 [2H] • 

Old reheating 

At the end of inflation, the inflaton fleld reaches the minimum of the potential 
and starts to oscillate coherently about this minimum. This is due to the friction 
term proportional to the Hubble rate in the Friedman equation f ll.lSp . We will 
now sketch some properties of the reheating scenario, since a general background is 
enough to deal with the content of this thesis. Let us consider a simple quadratic 
potential 



The inflaton is then described by sinusoidal oscillations, with decreasing amplitude 




(1.37) 



that follows from the approximated potential. 



1.3.2 Reheating 




(1.38) 



A{t), 



ip = A{t) sinmt , 



A{t) 



ITLp 



(1.39) 
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which give a decreasing energy density that is modified by the infiaton decay, 

p=^^^ + V{^) ^^m^A\t) oca-"^ . (1.40) 

In the old reheating model, there is a single infiaton field ip which is coupled to a 
scalar field x s-iid to a fermion ip, 

Ant = -o'V'X^ - h(fxl!^ , (1.41) 

with the coupling constants a and h. Here we have assumed that and are 
negligible with respect to the infiation mass m for simplicity. 

Quantum corrections are also taken into account, through the decay processes 
of ip. Let us call F = r((/9 — )■ xx) + r(v5 — )■ 'ipip) the total decay rate of the infiaton. 
The partial decay rates into scalar and fermion pairs result as [26] . 

n^^XX) = ^. r(^^^^) = ^, (1.42) 

which remain valid as long as F <^ m, namely cr^ ^ im? and <^ 1. By adding a 
phenomenological decay term Tp to the evolution equation, 

(^ + 3H0 + r0 + mV = 0, (1.43) 

we now obtain the solution 

A(t) = ^p-e-^*/^ (1.44) 



For small couplings a and /i, F < 3H at the beginning of infiation. Thus particle 
production may become important, in comparison to the total energy density, under 
certain conditions. This happens when 3H F, since the Hubble rate decreases as 
1/t. An estimation of the energy density of the Universe can be made by setting 
r2 = (3H)2 = 247rp/m2,: 

p = -. (1.45) 

By assuming that this is all transferred to light particles, which are instantaneously 
thermalized at temperature T^, we obtain 

«^^r^^ (1.46) 

^ 30 247r ^ ' 

is the effective number of degrees of freedom at T = Tr, and it is larger than 100. 
The estimation of the reheating temperature holds as 



T^<0.1v/fW, (1.47) 
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which by virtue of the relation F <^ m <, 10 ^mp can be rewritten as 

Tr<10-W- (1-48) 

This is a much smaller value than the temperature at the GUT scale Tqut ~ 
10~^mp ^ 10^® GeV. This means that the old reheating scenario does not provide 
a good GUT scale baryogenesis. 

Preheating 

We have seen that the old reheating model cannot describe correctly the phe- 
nomenology of particle creation in the early Universe. The theory of preheating 
[291 [30], on the other hand, is free of this problem. It postulates that before the per- 
turbative decay, the inflaton might have started to decay by means of a parametric 
resonance. This is a much more explosive process, which is non-perturbative. 

Preheating has a very complex dynamics, that is still not completely under- 
stood. Since it might affect significantly both Dark Matter production and the 
density porturbations, it has been studied extensively in recent years [311 [32l [33] . 
By including also this mechanism, the reheating process consists of three different 
stages: i) Preheating, where particles are produced non perturbatively by parametric 
resonance, ii) Perturbative decay of the inflaton, iii) Thermalization of the particles. 

We now consider this scenario in the hybrid inflation model, following Ref. 
[23] . The potential is given by 

r(^,x) = ^mV + ^^7W, (1.49) 

and we assume that spacetime and inflaton ip give a classical background, upon 
which the scalar field x is quantum. By expanding x in plane waves 

X = I (a.X.(t)e-^^-^ + 4x^(t)e^k-^) d'k , (1.50) 

and using the Friedman- Robertson- Walker metric, for each Fourier component Xk{t) 
the equation of motion is the following , 

Xk + SRxk + + Xk = 0. (1.51) 

By rescaling Xk and thus introducing the scalar field X^. = a^^'^Xk, the above equation 
is then recasted as 

Xk + ulXk = 0, (1.52) 
where the frequency of the mode k is 

+ (1.53) 
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During reheating, the term in brackets in (11.531) becomes neghgible, so it can be 
safely ignored. We then obtain from (11.521) the Mathieu equation 

^ + (Afc-2gcos2;2)Xfc = 0, (1.54) 
where z = mt. The amphtude of the resonance is 

A = 2g+4^, (1.55) 

with 

It is clear that the strength of the resonance in (11.541) depends on Ak and q. Regions 
of stability and instability are determined by the ratio Ak/q. In the unstable region 
there is production of Xk oc exp(/ifcz) and of particles with momentum k (/i^ is the 
Floquet index). For q smaller than unity, only a few modes grow and we have a 
narrow resonance. On the contrary, if q is much larger than the unity, resonance 
occurs for a broad range of the momentum fc-space. Since the growth rate of the 
produced particles is proportional to g, this is a much more efficient mechanism than 
the narrow resonance. It is called "broad resonance" [29] . 

From the above we see that the initial amplitude of the inflaton and the cou- 
pling g are fundamental to determine whether the resonance will be broad or narrow. 
Since the inflaton mass should be m ~ 10~^mp to satisfy the COBE normalizations, 
q is large ii g ^ lO"'^ with A{ti) 0.2mp. Then it can be shown that the broadest 
resonance is given by A^ = 2q. Eq. (ll.55p also implies that particles with low mo- 
menta are mostly produced. However, particles with high momenta can be created 
if g and Ak are large. This follows from the expression of the maximum comoving 
momentum |27] . 

t^^, (1.57) 

that is obtained by using the nonadiabaticity condition duk/dt ^ ul- For initially 
large q, there is a stochastic resonance of each mode [27]. In fact, the frequency cok 
decreases by cosmic expansion and changes within each oscillation of the inflaton. 
Therefore there is no correlation between the phases of the fields x ^i^nd if. This is 
important. In the first stage of preheating, the x fields cross many instability bands, 
without spending enough time on each band to provide with efficient resonance. This 
is in contrast to what happens in the case of Minkowski spacetime. However, as the 
cosmic expansion slows down, q becomes smaller and the fields stay in each resonance 
band for a longer time. Particle production ends when the variables decrease below 
the lower boundary of the first resonance band, by effect of the expansion of the 
Universe. 
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Interestingly, the resonance can be terminated by one more mechanism. There 
is a backreaction effect of the x particles which are produced, that modifies the 
equation of motion of the inflaton as follows, 

ip + 3Hip+ (m^ + g'{x')W = , (1-58) 
with the following expectation value of x^, 

(X^)^^ I k^\xk?dk. (1.59) 

An initial value q ^ 3000, which corresponds to g ^ 3.0 x lO"'', gives a growth 
of the variance (x^) of the order m? /g'^. In this case the backreaction becomes 
important p3j. This affects the oscillations of the inflaton field, which become 
incoherent, and the resonance is stopped. 



1.4 The matter content of the Universe 

In this chapter we have discussed models which aim to provide for a theoretical 
explanation of cosmological observations. It has been shown that a conservative way 
to explain the observed acceleration in the expansion of the Universe is assuming a 
certain form of the energy momentum tensor T^^ in the Einstein equations fll.2p . 

This corresponds to introducing exotic types of matter, such as Dark Matter 
and dark energy. In this final section we review the amount of each species as a 
fraction of the overall matter content, as it is obtained in agreement with recent 
experimental data. 

As we have seen in Section II. 2[ a key quantity to consider is the density pa- 
rameter VL. It is defined as the ratio of the actual (observed) density p to the critical 
density pc of the Friedman Universe. The critical density is the watershed between 
an expanding and a contracting Universe, since it corresponds to a static, fiat cos- 
mology with spatial curvature = 0. With this assumption, the first Friedman 
equation (II. 3p gives for the critical density today 

Pc(to) = ^ = ^^f^ = 1.88 X 10-^/.^^ , (1.60) 
87rG 87r cm'^ 

where G is the Newton's gravitational constant, Hq is the actual value of the Hubble 
parameter and the Planck mass is here expressed in grams as mp = 2.176 x 10~^g. 
The dimensionless parameter h is defined in term of the Hubble scale as 



H = 100/ikm/sec/Mps, 



(1.61) 



20 



CHAPTER 1. INFLATIONARY COSMOLOGY 



and its value has been recently measured as [20] 

/i = 0.73 ± 0.005 . (1.62) 

The density parameter of the species i is accordingly defined in function of the 
critical density, 

n - Pi SttGp^ 

^' = 7r^- ^'-''^ 

It is interesting to note that the above Pc(^o) corresponds roughly to 10 protons 
per m^. But actually the dominant matter is not baryonic, thus there are on an 
average only 0.5 protons per cubic meter! A general expression for Q where the 
density parameter equals exactly unity is commonly used. According to the ACDM 
model, the important components of Q are due to baryons (i.e. ordinary matter). 
Cold Dark Matter and dark energy. The WMAP satellite has measured the spatial 
geometry of the Universe to be nearly fiat, thus we have really k = and the total 
energy density is (almost) equal to the critical one: 

ntot = + ^DM + ■ (1.64) 

The total cosmological energy density is close to be critical, Qtot = 1 ± 0.02. This 
value has been obtained experimentally from the position of the first peak of the 
angular spectrum of the CMB radiation and the large scale structure (LSS) of the 
Universe. Let us now consider the different matter species separately. 

The ordinary, or baryonic matter, makes a very small contribution, Qb = 
0.044 ± 0.004, as found e.g. from the eights of the peaks in angular fluctuations of 
the CMB and from the production of light elements in the BBN. 

Dark matter makes instead a more relevant contribution, Qdm = 0.22 ± 0.04. 
This value is mostly made of two components: particles which are relativistic (called 
Hot Dark Matter), and those which constitute an almost inert relic in the intergalac- 
tic spaces, called Cold Dark Matter (CDM). There is also a third kind of particles, 
called Warm Dark Matter (WDM). WDM candidates usually are non-annihilating 
but rather weakly-interacting particles. They barely escape cosmological constraints 
like the free-streaming and self-damping bounds. In this case, structure formation 
occurs bottom-up from above their free-streaming scale, and top-down below their 
free streaming scale. They can be associated with any kind of particles, as long as 
they are just at the limit of the region allowed by self-damping and free-streaming. 

Dark matter is mysterious, since it interacts only gravitationally and as such 
it cannot be detected directly. However, its effects on galactic rotation curves, 
gravitational lensing, equilibrium of hot gas in rich galactic clusters, cluster evolution 



^We shall consider free-streaming for gravitinos in Chapter S) 



1.4. THE MATTER CONTENT OF THE UNIVERSE 



21 



and LSS provide with the actual value of ^dm- In this thesis, in particular, we focus 
on the possibility of producing Dark Matter gravitinos in the early Universe. 

Finally, the most mysterious and controversial component of the matter con- 
tent: dark energy. Nobody has a clear idea about what it really is, either a purely 
repulsive form of energy or some peculiar scalar particles (moduli, quintessence 
fields...). It drives the accelerated cosmological expansion and is uniformly dis- 
tributed in the Universe. The dimming of supernovae with high redshift, LSS, the 
CMB spectrum and the age of the Universe give the value fl fa 0.74. Since this 
is overwhelmingly the dominant component of the total matter content, we say 
that the actual Universe is dominated by a cosmological constant, in contrast with 
previous periods of radiation and matter domination. 
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Chapter 2 

About supergravity and gravitino 
cosmology 

In this chapter we consider supergravity (SUGRA), namely the theory that corre- 
sponds to local supersymmetry. The gauge field is the gravitino, which obtains a 
mass through the super-Higgs mechanism [SU [351 EE] . Our discussion is based on 
the review by Lyth and Riotto ^14j and on Takeo Moroi's thesis [37], to which we 
refer for deeper analyses. 

2.1 Basics of supersymmetry and supergravity 

The Standard Model of Particle Physics (SM) is an established and successful frame- 
work, which has been tested experimentally with great accuracy. Nevertheless, it 
is common conviction that it has to be somehow extended to more fundamental 
theories. There are not really inconsistencies within the SM, however the theory 
does not seem completely natural (remember the famous " hierarchy problem" ) . 

To this aim, supersymmetry (SUSY) is certainly the most promising possibil- 
ity. In very simple words, it interchanges SM (or "ordinary") particles with new 
particles, which have higher masses and are called "superpartners". SUSY provides 
with a natural solution of the hierarchy problem and it gives a remarkably good uni- 
fication of the gauge constants. Moreover, via conservation of R parity, it introduces 
automatically a candidate for Dark Matter. 

We still lack experimental evidence of the superpartners, thus SUSY must be 
broken to account for the mass difference between the two sets of particles. In the 
case of global supersymmetry (the transformations are not space-time dependent), 
the SUSY breaking cannot be phenomenologically satisfactory. For example, there 
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is no experimental candidate for the massless spin-1/2 particle called goldstino, that 
is created by the global SUSY breaking. 

As a solution to this drawback, one might consider local transformations. The 
gauge particle is the spin-3/2 gravitino, which after the breaking of supersymme- 
try absorbs the goldstino modes and becomes massive. Since SUSY requires the 
gravitino to be coupled to a spin-2 ordinary particle, we find this in the graviton. 
Gravity is then automatically embedded, and supergravity (SUGRA) is the theory 
which arises accordingly. This provides with new interesting insights, from both 
theoretical and experimental perspectives. 



2.1.1 The super symmetry algebra 

SUSY theories of type N=l (see for instance [35], EHl [39]), which have only one 
fermionic generator, take part to the lowest class of supersymmetric models. Yet 
they are relevant for inflation. In the following, we describe their basics with the 
conventions of Wess and Bagger [36] • 

Let us consider a supersymmetry algebra given by 

{Q^,Q^} = 2a^^.P,, (2.1) 

where Qa and are the supersymmetric generators and the bars stand for con- 
jugation. The two-component Weyl spinors are labelled by a and /3, which run 
from 1 to 2. The quantities with dotted and those with undotted indices transform 
respectively under the (0, |) representation and the (|,0) conjugate representation 
of the Lorentz group, cr'^ is a matrix four vector, cr^ = (— 1,(t) and is the 
four-momentum operator. 
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The supersymmetry algebra that contains fields of spin less than or equal to 
one have two irreducible representations, the chiral and vector superfields. More 
precisely, chiral fields contain a Weyl spinor and a complex scalar, vector fields 
contain a Weyl spinor and a massless vector. A chiral superfield can be expanded 
in terms of the Grassmann variable 9 in supcrspacc, 

0(a;, 6) = 0(x) + V2eilj{x) + e^F{x) , (2.2) 

where x is a point in spacetime, is the complex scalar, ip the fermion and 

F is an auxiliary field. Following a common habit, here we maintain the same 

symbol to represent a superfield and its scalar component. Under a supersymmetry 
transformation with anticommuting parameter e, the component fields transform as 

d(f) = V2eijj, (2.3) 
5^jJ = V2eF + V2ia''ed^(l), (2.4) 
5F = -V2id^ijja''e . (2.5) 

On the other hand, the vector superfields in the Wess-Zumino gauge, for the simplest 
case of an abelian group U(l), can be written as 

V = -Oaf'eA^ + ieWx - iPex + \qWd , (2.6) 

where A is the complex conjugate of the two-component Weyl spinor A. Here A^ 
is the gauge field, A is the gaugino and D is an auxiliary field. The analog of the 
gauge invariant field strength is a chiral field, 

W^a = -^Aa + - \{a^-a''e)o.F^. + ^V^^^A^ , (2.7) 

where F^i, — d^Ai, — di,A^ and — (—1, — <?)■ Under the supersymmetry transfor- 
mations, the gaugino A transforms as 

bX = ieD + ea^'d'^F^^ . (2.8) 

Global supersymmetry corresponds to invariancc with respect to SUSY transforma- 
tions with parameter e independent of spacetime coordinates. In local supersym- 
metry (i.e. in supergravity) , e is instead x^-dependent, thus we write e{x). In the 
latter case, one has to introduce one more supermultiplet that contains the graviton 
and the gravitino. 

Global SUSY can be seen, with a rough estimate, as the limit of supergravity 
when the Planck mass goes to infinity. However, this is a good approximation only 
if the vacuum expectation values (vevs) of all the relevant scalars, which have not 
been integrated out, are much less than mp. 
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Moreover, the overall picture is not completely clear. In the true vacuum, 
global SUSY would predict a large positive potential V, while the observed value 
is almost zero. Siipergravity instead predicts V ^ 0, through subtraction of an un- 
known large value. This generates the cosmological constant problem. The second 
exception deals with the limit mp — > oo, that does not make any sense during infla- 
tion. The Planck mass is indeed a fundamental parameter during that epoch. The 
naive limit of an infinite Planck mass should be replaced by less immediate mech- 
anisms. This generates a problem for model building, since supergravity theories 
imply in general an inflaton potential that is not sufficiently flat for inflation. 

By keeping this in mind, we now discuss the Lagrangian of a globally super- 
symmetric theory. The fundamental concepts of superpotential, scalar potential and 
F and D terms will be now introduced. 

We begin with the most general renormalizable Lagrangian written in super- 
space, 

>C = ^ y ci^^ 0i,e^0n + ^ / + J W{<Pn) + h.c. (2.9) 

with the superpotential W{4>n{x, 0)) and with the matrix M = T"V". A vector su- 
perfield V" corresponds to each generator of the gauge group G, in the represen- 
tation that is determined by the In the adjoint representation, Tr(T'"T'^) = kS'^^. 
Consider now for simplicity the gauge theory of a single U(l), with coupling constant 
g and charges Qn- The covariant derivative of a chiral superfield, 

g 

D^(t)n = df^cpn - i-A^qn(t)n , (2-10) 
allows to rewrite the Lagrangian in terms of the component fields as follows: 



jC = 



n 



^/2 



1 d'^W 

2 d(f)nd(f)r 



dW 

d(t)n 



-t- c.c. 



(2.11) 



Clearly, D and F„ do not propagate, since they do not have kinetic terms in the 
above equation. This implies that they are auxiliary fields and their equations of 



2.1. BASICS OF SUPERSYMMETRY AND SUPERGRAVITY 



27 



motion are just simple constraints, 

D -- 

The above equations f l2.12p and fl2.13p can then be used to extract the scalar po- 
tential from the Lagrangian Eg. ( 12 .11]) . and to separate it into the F and D term as 
V = Vp + Vn, where 

Vf ^ (2.14) 

n 

Vd = \d\ (2.15) 

In any renormalizable theory, W is at most cubic in the fields (thus the potential 
is quartic), though in several cases higher order terms are admitted. We will see 
below that the above splitting is extremely important in any generic SUSY theory, 
including super gravity. 

Terms like {l/2)m(f)f in W are usually forbidden in theories which arise from the 
Lagrangian dMUlll. This means that the masses in the observable sector come only 
from SUSY breaking and from the spontaneous breaking of an internal symmetry. 

The famous Fayet-Iliopoulos term [ID] 

-2^ j d^9V, (2.16) 

corresponds to adding a contribution — ^ to the D field. Eq.f l2.13p thus becomes 

^ = -|E^"l'^-l'-^ (2.17) 

n 

The D term is accordingly modified as 




(2.18) 



If ^ and the charges are redefined so that 




(2.19) 



"'^With the only exception of the /i term iiH^H^, that gives mass to the Higgs doublet. 



fdW\ 

' \Wn) 

9 



i:^qn\(l)n\^ ■ 



(2.12) 
(2.13) 
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The additional term ^ is usually related to internal symmetries, and it can appear 
in the theory from the very beginning. For instance, weakly coupled string theories 
allow anomalous U(l) symmetries which provide with terms such as fl2.16p in the 
Lagrangian. 



2.1.2 Global supersymmetry breaking 

If supersymmetric models aim to be phenomenologically realistic, supersymmetry 
must be broken. In fact, we did not observe so far any superpartner of the parti- 
cles we could detect at colliders. This means that the ordinary particles and their 
superpartners are not degenerate in mass. The breaking of global supersymmetry 
may occur either explicitly (i.e., by adding soft breaking terms in the Lagrangian) 
or spontaneously. 



Soft SUSY breaking 

Spontaneous breaking of supersymmetry can be done in several ways, and is still 
object of discussion. The main idea behind the explicit, or "soft" breaking is instead 
very simple: since in the effective theory at low energy (below 1 TeV) SUSY is broken 
explicitly, we just add to the effective Lagrangian density some SUSY-breaking 
terms [H]. We do not worry about their origin. However, the choice is somehow 
constrained, since the naturalness of the theory admits only certain soft breaking 
parameters. 

The couplings must be "soft", namely they have positive mass dimension to 
ensure a natural hierarchy between the electroweak and the Planck scales. These 
can typically be the scalar and gaugino masses or the cubic terms in the scalar 
field potential {A^^'^cpiCpjcpk, with couplings A^^^). Soft quartic terms vanish in a fiat 
direction of unbroken SUSY. In fact, phenomenology requires the couplings and the 
scalar and gaugino masses to be smaller than 1 TeV |12] . 

The squark and slepton masses come almost entirely from the supersymmetry 
breaking (with the exception of the stop). They should be 100 GeV, since they 
have never been observed experimentally. Such considerations give for the masses 
of the scalars m an expected range between 100 GeV and 1 TeV. The possible soft 



we obtain 




(2.20) 
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SUSY breaking terms in the Lagrangian are [H] 

Csoft = - Qm, A'^A'^ + ^A'^'Mj<f)k + ^M'^Mj + t'^?j + c.c. - (m^);0^*0, . 

(2.21) 

We recognise here the gaugino masses Ma for each gauge group, the scalar (squared) 
mass terms M*-^ and (m^)*, and the tadpole coupling f. The latter is associated 
with a gauge singlet, and as such it is absent from the MSSM. In Ref. [33] it is shown 
that the above soft breaking terms do not generate any quadratic divergences from 
quantum corrections to the scalar masses, to all orders in perturbation theory. 

The terms in Eq. (12.211) break supersymmetry because they contain only scalars 
and gauginos, not their respective superpartners. Thus they give masses to all the 
scalars and gauginos of the theory, even though the gauge bosons and the matter 
fermions are massless or very light. Gauge symmetry allows both Ma and (m^)*. 
The other couplings instead are allowed by gauge invariance only if there is a corre- 
sponding term in the superpotential. 

The soft breaking terms are put in the Lagrangian by hand. As already re- 
marked, we do not know how they are generated unless we invoke the spontaneous 
supersymmetry breaking. Generally speaking, the effective theory, which is valid 
only below the breaking scale, describes uniquely the observable sector, which con- 
tains the fields of the Standard Model and their supersymmetric partners. 

In the full theory the spontaneous SUSY breaking occurs in a hidden sector, 
where the fields do not obey to SM gauge interactions. The effective theory in the 
observable (or visible) sector is obtained when the hidden sector is integrated out. 
This happens when the scalar field potential, at fixed values of the fields which 
remain in the action, is minimised with respect to the hidden sector fields. This 
gives a well-defined theory if their masses are much larger than those of the fields 
that are not integrated oulU (or if the two sets are coupled only gravitationally), as 
it is discussed in [14] and in [35j. 

How the SUSY breaking is mediated from the hidden to the observable sector 
defines different classes of theories. If this happens through gravity, the theory 
is called "gravity-mediated". On the other hand, if the mediation occurs through 
gauge interactions, we talk about "gauge mediation". 

Fayet-Iliopoulos supersymmetry breaking: the D-term 

We now investigate the case in which supersymmetry is spontaneously broken. This 
means that there is still invariance (or covariance) of the Lagrangian under SUSY 

^In fact, in this case their motion about the minimum of the scalar potential is negligible. 
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transformations, but there is no invariance of the vacuum. More into detail, the 
generators Qa do not annihilate the vacuum. They produce instead a chiral fermion 
ipa or a gaugino A^. This happens if the vev for {Qa, V'/?} or {Qa, A/3} is non-zero. 
The transformation rules of the chiral fermion (12. 4 p and of the gaugino (12. 8p define 
these quantities. SUSY is spontaneously broken if at least one of the auxiliary 
fields Fn or D has a non-vanishing vev. In the true vacuum, the scale of global 
supersymmetry breaking is defined as 

^^ = J2\Fn\' + \D' = V. (2.22) 

n 

Spontaneous supersymmetry breaking can occur either at tree-level or dynamically 
(namely, due to quantum corrections). In any case, the mass spectrum is modified, 
since the symmetry between scalar and fermion masses is lost |4T1 144] . In the simplest 
case, D = and there is just one Fn. This is the O'Raifeartaigh model, that we will 
discuss in the next subsection. Here we focus on the alternative D ^ 0: the D term 
drives the breaking mechanism. This corresponds to the " Fayet-Iliopoulos model" . 

Consider a non-zero D-teim in the SUSY Lagrangian, which we write as [H] 

CFi = -kD. (2.23) 

The dimension of the constant k is [massY . Since this term is linear, we expect D 
to get a non vanishing vev. By recalling Eq. (I2.20p , we can write the relevant part of 
the scalar potential as 

V = kD-]^D^-gDY, qn\<Pn? . (2.24) 

n 

This gives the equation of motion 

D = k-gJ2Qn\(l>n\\ (2.25) 

n 

therefore if the scalar fields 0„ (which are charged under the U(l) group previously 
introduced) all have non-zero masses m„, the scalar potential becomes 

V = J2 l^^nlVnl' + lik-gJ2 Qn\^n\y ^ . (2.26) 

n n 

Since the minimum occurs for non vanishing D, there is spontaneous supersymmetry 
breaking. The scalars have squared masses (|m„p — gqnk), whilst for their fermion 
partners we have |m„p. 

It seems that the above mechanism provides spontaneous SUSY breaking at the 
tree level in a rather simple fashion. Anyway, there is a fundamental problem with 
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the Fayet-Iliopoulos term. Up to now we have considered only an Abehan group. For 
non-Abehan Yang Mills theories, a term like f l2.23p would break gauge invariance. 
It follows that only D-terms in U(l) can be used, which is bad for phenomenology. 
In the MSSM, in fact, the squarks and sleptons do not have mass terms in the 
superpotential, so they would get non-zero vevs in order to make the equation of 
motion (12.251) vanish. This would break colour symmetry, but not supersymmetry. 
Accordingly, a F-I term must be subdominant with respect to other mechanisms of 
SUSY breaking, if not totally absent. 

Therefore we can conclude this discussion by saying that, even though the 
D-teim breaking has not been completely ruled out nowadays, it does not seem 
to be efficient in giving the right masses to the MSSM particle spectrum (and in 
particular, to the gauginos). The F-term, or O'Raifeartaigh model, gives instead a 
more transparent and less problematic phenomenology. 



O'Raifeartaigh supersymmetry breaking: the F-term 

If only the F-term is involved, the corresponding models of tree-level SSSB are 
called of the O'Raifeartaigh type [15]. The SUSY breaking is due to a non-zero vev 
of the F-term. The superpotential is chosen so that the equations Fj = have 
no simultaneous solutions and the scalar potential V = |Fjp is positive at its 
minimum. More into details, let us consider the set of three chiral supermultiplets 
{S, X, Y) and the superpotential 

W = -fi^S + XSX^ + mXY , (2.27) 

which we will discuss further in Chapter HI The key element in W is the term 
—/Ji^S, which is required in renormalizable theories with tree-level SUSY breaking, 
if the superpotential contains only polynomial terms |l6]. In fact, since A and m 
can be assumed to be positive with no loss of generality, the term —fi'^S implies the 
equations 

Fs = ^ = -fi' + XX^ = 0, FY = ^ = mX = 0, (2.28) 

which are incompatible. The equation for X, 

dW 



dX 



2XSX + mY = 0, (2.29) 



is nothing more than a constraint on the fields. Therefore supersymmetry is broken, 
and the scalar potential is the following: 



V{S,X,Y) = 



(2.30) 
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The absolute minimum of is at X = F = 0, with S left undetermined. This 
means that at the minimum 

Fs = -^"^ = (2.31) 

Thus fi is the SUSY breaking scale at the tree level, and S' is a flat direction of the 
scalar potential, namely a subset in the scalar field space along which the potential 
vanishes (by rescaling). In the MSSM, around 300 such fields have been catalogued 
[47] . Flat directions are interesting also in the context of cosmology, since they 
might explain the origin of the inflaton [ilSj. 

Let us now consider the particle spectrum of the O'Raifeartaigh model. If V 
is expanded around S = 0, there are 6 real scalars with squared masses 

0, 0, m\ m\ - 2A/i^ + , (2.32) 

where the eigenvalues correspond to the complex scalar 05, which is the scalar 
component of the chiral superfield S: 

S{x, 6) = M^) + V29Mx) + O^Fsix) . (2.33) 

The fermionic sector contains three Weyl fermions with squared masses 

0, (2.34) 

thus the fermionic partner i/^s oi the massless scalar (^'5 is massless as well. It is 
anyway important to remark that ips exactly massless, since it is the goldstino. 
Instead, the scalar is massless because of the fiat direction. However, quantum 
corrections remove ("lift") the fiat direction, and generate the one- loop Coleman- 
Weinberg potential which provides the scalar 0s with a mass 



2 \2 2 



2x , , 1 / . 1 \ . fl + r 



(2.35) 



ln(l - r^) - 1 + - + - J In ^ ^ 

where r = 2A/i^/m^ (this mechanism will be considered explicitly in Chapter H]). We 
see from the above equations that the two sectors are not degenerate, as it follows 
from the spontaneous breaking of supersymmetry. It is straightforward to check 
that the general supertrace formula, written as 

'-ifml = 0, (2.36) 



where (— 1)"^ is 1 for bosons and —1 for fermions, is verified in this model. This is 
a general result of renormalizable theories at tree level. This simple model seems 
already satisfactory, but there is a free parameter, namely the SUSY breaking scale 
II, which needs to be fixed by hand. Crucially, fi must be much smaller than mp 
to be consistent with the phenomenology of the weak interaction. The scales of the 
soft terms are generated naturally only in theories with dynamical supersymmetry 
breaking, where they appear by virtue of dimensional transmutation [IHl SB]- Such 
a mechanism will not be addressed in this work. 
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2.1.3 The general supergravity Lagrangian 

So far we have considered only local supersymmetry. The SUSY algebra and soft, D- 
term and F-term breaking have been defined and discussed. However, the gravitino, 
which is the main topic of this thesis, does not exist until one introduces supergravity. 
This is what we are going to do in this section. 

As it is well known, gravity is a non-renormalizable field theory. The same 
holds for supergravity, which by construction is induced by gravity. As such, it 
is regarded as an effective theory, that is valid only below a certain ultraviolet 
cutoff Auv In the case at hand, this scale is the (reduced) Planck mass Mp = 
2.44 X lO^s GeV. 

The supersymmetry transformations, which are defined by Eqs.f l2.3l) . (12. 4 p and 
(12. 5p . still remain valid in supergravity. The difference is found in the Lagrangian, 
since in addition to the superpotential W, we now define two more functions: the 
Kahler potential K and the gauge kinetic function /. W and / are holomorphic 
in the complex scalar fields, but K is not. In this context, the only physically 
meaningful object is the so-called Kahler function 

We stress that within supergravity, G and / are completely arbitrary. The choices 
which are made in the literature are indeed based on purely phenomenological 
grounds. 

Before writing down the Lagrangian, let us consider the above three fundamen- 
tal functions, namely the Kahler potential K, the superpotential W and the gauge 
kinetic function /. The Kahler potential is a real-valued function K{(f), cf)*) in the 
chiral scalar fields, with mass dimension two. Since it is not a holomorphic function, 
it is not strongly constrained by symmetries. It determines the kinetic terms of the 
scalar fields as follows, 

A.n= (5^0")^7.*i(9'^0"), (2.38) 
where we have introduced the Kahler metric (7, whose components are defined as 

= » - '2-3^' 

along with its inverse g"^^ , namely g^-^* gj*k = 5^- The Kahler connection is given by 

r& = /'-^- (2-40) 
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In the case of canonical Kahler potential, the scalar fields are canonically normalized 
at the origin. This corresponds to gij* = 6ij*, and K can be written in the following 
generic form, 

K = Kh{h, h*) + aij{h, h*)(f)'(f)*' + [Zij{h, h*)(f)'(f)^ + h.c] . (2.41) 

The hidden part of K that is independent of observable fields has been labelled 
with Kfi. The above definition, along with the Kahler metric, defines an interesting 
topology which is based on the concept of Kahler manifold. The above metric and 
connection define indeed the covariant derivatives of the superpotential, 

DiW = Wi + MfKiW , (2.42) 
DiDjW = Wij + Mp\KijW + KiD.W + K^DiW) - T^^DkW + 0{Mf) , 

(2.43) 

with the compact notation Ki = diK, Wi = diW, Kij = didjK, and Wij = didjW. 

The next object we consider is the superpotential W. As already discussed, it 
is an analytic function in the chiral scalar fields, and it has mass dimension three: 

W = Wh{h) + i/i.,(/i)0V'' + lmjk{hW<j>><P'' . (2.44) 

In contrast with K, the superpotential is holomorphic, thus it must be consistent 
with gauge invariance. Here we have left the superpotential of the hidden sector 
Wh{h) implicit. The simplest option for ly is a direct sum of the superpotentials of 
the hidden and observable sectors, 

W = Wh{h) + WoU<P), (2.45) 

even though in theories with dynamical SUSY breaking, a mixing of the two sectors 
through a hidden field fiij{h) is generally required [46] . 

The gauge kinetic function / determines the kinetic terms of the gauge fields 
and of the corresponding gauginos. Since it is holomorphic, internal symmetries 
constrain its form exactly as it happens with the superpotential. / is dimensionless 
in the scalar fields 0, and it multiplies the kinetic term for the vector supermultiplet. 
This means that the derivatives of fab{4>) with respect to the scalars 0*, 

dja, - ^ , (2.46) 

have negative mass dimension. This is relevant in theories of spontaneous SUSY 
breaking, as in certain models difab = C(Mp^). Accordingly, terms in the super- 
gravity Lagrangian fl2.66p which are proportional to {Mp^ difab) are C(Mp^), thus 
sub dominant. 
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The general supergravity Lagrangian has a number of symmetries. It is invari- 
ant under the following transformations, where F{(f)) is an arbitrary holomorphic 
function, 

Ki<P, 0*) ^ 0*) + + , (2.47) 

W ^ e-^^^^/^^'^W , (2.48) 

provided the following F-dependent Weyl rotations of the spinor fields: 

^ ei^-^/^'^^^^xl, (2.49) 

A« ^ g-flmF/Afl^s^a^ (2.50) 

^ e-i^-^/*^^^^^^. (2.51) 

Eq. fl2.47p leaves the metric Qij* and the Kahler connection T^j unchanged. The 
isometrics of the Kahler manifold are described by the real Killing potentials, which 
we call Da{(f), 4>*). These are found by solving the corresponding Killing equation in 
the Killing vectors, 

Xi = -zg'^*d,,D,, (2.52) 
=tg'^*d,Da, (2.53) 

which obey the field transformation 



^ (f)" = (j)' + X''\(f))e , (2.54) 
^ = (f)i* + X*''\(j)*)e, (2.55) 



with the infinitesimal parameter e. The Killing equations 

ViX; + V,Xf = , (2.56) 
V.Xf +V,-.Xf = 0, (2.57) 

follow from the explicit form of the Lie derivative Cx of the (Kahler) metric Qij*: 

~ Q d ^ d J, 

= VjX| + VjXf = 0. (2.58) 

In the above equations, X^ = Qi-jX""^ and the index i represents both i and i* . The 
covariant derivative V of a generic "vector" Ai has the generic form 

V.A, = Aa^. _ T\^A, . (2.59) 
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Eg. (12.561) is identically satisfied, while fl2.57p allows to express the Killing vectors 
and X*°-'^{(j)*) as derivatives of the Killing potential D"^: 

X-{^) = -^g^^'^D\ (2.60) 
X*''\(j)*) = tg^rA^D" . (2.61) 

This is interesting, since the above can be used to define an analytic function F"" of 
the scalar fields as follows. 

This object appears in the covariant derivatives of the fields, which we will address 
later. For example, with a minimal Kahler potential i^*"*" = <;/)*(;/)** and the genera- 
tors of a gauged Lie group T^^-, the Killing vectors and the Kilhng potential are the 
following: 

X^\<1>) = (2.63) 
X*^\(f)*) = i(j)*^T^i, (2.64) 

^ (f)*%^.(f)^ . (2.65) 

By using the Kahler potential, the kinetic function / and the superpotential W, 
the general form of the supergravity Lagrangian can be now obtained. As a formal 
remark, we point out that the following is an all-order on-shell Lagrangian. This 
means that the auxiliary field of the SUGRA multiplet is eliminated via the field 
equations of the vielbein. In fact, the kinematical null-torsion constraint = 0, 
that is used in this standard approach, allows to express the spin connection u'^^ in 
function of the vielbein and of the gravitino ip^. Therefore the spin connection is 
no longer an independent degree of freedom [SU HHl EO] • 

In the following equation we impose Mp = 1, = Re/ and = Imf and 
(a, b,c, ...) are indices of the adjoint representation of the gauge group with coupling 
constant g and structure constant f°'^'^. The indices are raised and lowered with 
and its inverse. The Lagrangian is very complicated, and it contains the scalar fields 
0, chiral (matter) fermions gauge bosons AJJ and gauginos A, the vielbein and 
the gravitino ^pf, [32 [Ml E?]. 
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--^/2egdJa,D'^x')^'+\V2egd,*f:,D'^r>^ 

+^egDaiP^a^'~X'' - ^^egDai^^^a^X" 

-^V2egij*D,(l)*^X'c7''a''il:^ - ^V2egij*D,(l)'x'^''(^''^,. 

--e {ij.aVXa + i^.aVXa) {f^, + F^^ 

+^eg,j. {ie^'''"^ij,a,^, + ^.a'^r) x'^cX' 

-^e {gij*gki- - 2Rij,ki*) x'x^X^X^ 

+ [2g,^.f^, + f'^'^'^-'d^h^d^.taa] x' ^^x'~X^^,X' 

+^eV,a,/„a¥A"A^ + ^eV,.a,-./:,x¥A"A^ 

-^e^'^*a,/„fea,-./:,A"A^A^A'' + ^eA.a'^A^Aba^A'' 
+ ]y2ed,h, [x^a^^'^X^^l^^a^X' - ^V^^a'^x'A^A^^ 

+ '-V2ed,*f:, (r^^'-Xy.aA' - ^V^a'^x'A^A^^ 

+^V^ee^/2 {D.Wx'a^^l^, + D,.W*x'^^^,] 
—ee""^ {D,D,WxY + D,,D,,W*x'x'} 

+^ee^/^^^^* {Dj.W*dJabX^X'' + AW^9,*/:,A«A^} 

-ee^ [g'^\DiW){Dj.W*) - 3W*W] . (2.66) 
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The covariant derivatives are defined as 

^ d,ct>'-gAlXl, (2.67) 

{K,D^<P^ - K,,D^ct>*^) X' - '-gAllmFax' , (2.68) 
+ ^ {K,D^4>^ - K,,D^<P*^) \- + '-gA'^lmF,r , (2.69) 

+ ^ {K,D,^ - K,,D,^*^) + ^(^^ImF.V. • (2.70) 
is the field strength (or curl) of the gauge field A^, and is defined as 

F% (^'^"^'^^^ + + ^'^^'^^^ + ^-^'^^") ■ (2.71) 

The supergravity Lagrangian (12.661) is invariant under local supersymmetry 
transformations. These are parametrised by an anticommuting Majorana spinor ^ 
of mass dimension —1/2, that will be written here as ^(a;), with abuse of notation. 
For the bosons (the vielbein e^, the scalars (ff and the gauge fields A^), we have 

5e/ = -i{ia'^i^, + icy'^ij,) , (2.72) 

6^' = V2^x\ (2.73) 
= ^(ecr^r + e>^A'^) . (2.74) 

For the matter fermions x\ the transformation rule is the following, 

Sx' = %^o^ {d,<^' - ^^mX^) - ^)M'X^ + \ {KM' - K,^'') X' - 

-V2F^i + A^ , (2.75) 

and for the gauginos A" one finds 

+ ^^f''^'''^-%fMX'\' - ^e/''^"''^-'5.-/(;e)X^A^ • (2.76) 
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Finally, for the gravitino i/j^^ we obtain 

_1 _ Kj*6(l)'') + te'^/^Wa^^ . (2.77) 

The auxiliary field F\ that plays a crucial role in global SUSY, here appears by 
virtue of the following definition: 

F' = e^/^g'^*Dj*W* . (2.78) 

Let us now focus on the particle content of the Lagrangian, that consists of 
three distinct sectors. Matter fermions are described in terms of left-handed four- 
spinors 

Xl = (^0^) , (2.79) 

with the two-component Weyl spinor Xa- The scalar superpartners are denoted as 
and the multiplet is in the fundamental representation of the gauge group. The 
gauge bosons and the gauginos 

A" = t ("[Is) , (2.80) 

namely Majorana fields which are their superpartners, constitute the gauge super- 
multiplet, in the adjoint representation of the gauge group Q. The indices a,b,... 
run accordingly. Since SUGRA is a Yang- Mills theor}0, F^^^ are the associated field 
strengths of the A^. Moreover, the auxiliary fields Da are the generalizations of the 
D terms in the vector supermultiplets in global supersymmetry. 

The gravity multiplet is a bit different. The graviton (which we indicate as 
the vielbein ej^), enters the SUGRA Lagrangian via the determinant of the vielbein 
e = dete™. It also appears in the curvature (Ricci) scalar R. By recalling the 
fundamental relation 

9,u = Vmne'^el, (2.81) 

we distinguish between the flat spacetime indices (m, n, . . . ) and the Lorentz indices 
(/i, z/, . . . ). We will consider the gravitino into details in the next section. It is a 
spin-3/2 fleld, which is written in terms of the Majorana vector-spinor, 

= ^ (""^j) • (2.82) 

The gravitino is massless in the limit of unbroken local SUSY. As we will see in the 
next section, the SUSY breaking provides it with a (f ef-dependent) mass, through 
the super-Higgs mechanism (see for instance Ref. |37]). 

•^N^l supergravity is indeed the gauge theory of the Poincare supergroup IS0(4,1). 
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2.1.4 Spontaneous breaking of local SUSY 

In contrast to global SUSY, supergravity can be broken only spontaneously, not 
explicitly. However, the condition for the breaking is still a non-vanishing vev of 
at least one of the auxiliary fields D and F^. The similarity with the spontaneous 
gauge symmetry breaking of the Standard Model is however less evident in this case. 
The vevs of the D and can indeed get contributions not only from scalar fields, 
but also from fermion condensation (for instance, from gaugino condensation). 

The scalar potential at the tree level is given by 

V = Vd + Vf, (2.83) 

similarly to the case of global SUSY. However, here Vp is not only determined by 
the F-terms, because the Kahler potential makes the difference. 



|2 



(2.84) 



The ultraviolet cutoff k\jy can be now recalled and combined with the reduced 
Planck mass Mp = 2.43 X 10^^ GeV. Taking Mp to infinity and keeping the cutoff 
fixed, we have 

Vp = Wng"''''iWmr. (2.85) 

Namely, we have obtained non-renormalizable global supersymmetry. Renormaliz- 
able global SUSY is recovered if also Auv goes to infinity. If instead Auv oo 
and Mp stays fixed, the theory corresponds to minimal supergravity, which contains 
only the graviton and the gravitino '0^. In such a model, the total Lagrangian 
is the sum of the Hilbert-Einstein and Rarita-Schwinger Lagrangians (the latter is 
basically the kinetic term of the gravitino field): 

C = Che + Crs = -^eR + ee^^'^P^i^^a^Dpi;, . (2.86) 

Cps ensures invariance of the action under local SUSY transformations, and it is 
for this reason that the gravitino is the gauge field of local supersymmetry. Usually 
we consider Ajjv = Mp, as in the Lagrangian (I2.66p . Global supersymmetry is 
recovered when the cutoff (namely, the reduced Plank mass) goes to infinity. 

We now study into detail the spontaneous breaking of supergravity, following 
[35] and [37]. As already discussed in the previous sections, supersymmetry is broken 
by vevs acquired by fields in the hidden sector, which are gauge singlets. This 
generates soft breaking terms of order 



(2.87) 
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which vanish in the hmit of exact supersymmetry. That is, when the vev of the 
auxihary field approaches zero, (F) — )• 0. 

It is very important to remark that the different ways of breaking SUSY 
(through gravity, gauge interactions,...) do not change the structure of these pa- 
rameters. What is different in each model is their size, which determines the mass 
spectrum of the various particle multiplets. Thus the Lagrangian (12.661) is completely 
general. Accordingly, this occurs also for the interaction terms and the Feynman 
rules which are used throughout this thesis. 

After that global SUSY is spontaneously broken, a massless Nanbu-Goldstone 
fermion called the goldstino is generated. In local supersymmetry, this component 
is absorbed by the gravitino, which obtains a mass (super-Higgs mechanism). The 
gravitino mass is therefore directly related to the SUSY breaking scale. Let us 
refer once again to the general supergravity Lagrangian Eq. (l2.66p . In the following, 
we will consider only the minimal (or canonical) Kahler potential and the kinetic 
function at the lowest order, namely 

ir = ^0>« + ..., f(ab)=Sab + ... (2.88) 

i 

which is a good enough approximation for the purposes of this thesis. The dots 
in the above equation label higher order terms in the inverse reduced Planck mass 
Mp, which describe subdominant interactions. We remind that K and / generate 
respectively the kinetic terms of the chiral and of the vector multiplets. The choice 
(12.881) of a canonical Kahler potential means that gtj* = 6ij* + . . . , namely it defines a 
flat Kahler manifold. Extracting from Eq. (12.661) the quadratic fermion terms without 
derivatives [37] , 

+ i (G,.^.. + x^x^j + V2g {-tDtx'X^ + tDl^x''^'^) , (2.89) 

we notice mixing terms of the gravitino with either the matter fermion or the 
gaugino A*^, provided that the vevs of e^/'^DiW or D'^ are nonzero. It is possible to 
eliminate such terms by a shift of the gravitino field, 

^ V^M + , (2.90) 
where the spinor 77 is defined such that 

r?^^G..x^ + ie-^/2^Z}'^A". (2.91) 
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The new form of the Lagrangian is 

(2.92) 

By restoring the exphcit dependence on the Planck mass, we accordingly find the 
gravitino mass as 

m3/2 = = -^\W\. (2.93) 

It can be shown that the spin 1/2 fermion is a massless eigenstate of the fermion 
mass matrix in the basis (x*, A"), thus it is the goldstino. We have thus verified that 
this particle is absorbed (or "eaten") by the gravitino field ip^. 

We conclude this section with a discussion on the supertrace formula in super- 
gravity, that was here introduced in the context of global SUSY, Eg. (12. 36 p : 

3/2 

StrM^ = ^(-1) . (2.94) 

j=o 

Here each Mj is the squared mass matrix of the corresponding particle multiplet, 
namely the set of scalar, vector and spinor fields. The mass matrix of the scalar 
fields can be obtained from the second derivatives of the scalar potential V , in the 
case of a vanishing cosmological constant. The trace is then 

trM^ = (_n^ -1) _ie-^(/D«D«)V2e^(G,,G,*,.+n^) 

+2g\DtD%+D'^Dl,) , (2.95) 

where = Y^- gu* is the number of chiral multiplets [3^ . The mass terms for the 
vector bosons are derived from the covariant derivatives of the scalar fields. The 
trace of the mass-squared matrix is 

trM? = Gg'^T^^T^,^' = Gg'DtD^, , (2.96) 

where the Killing potential D"^ for the minimal Kahler potential is obtained from 
Eq. (l2.65|) . Finally, the super-Higgs mechanism provides with the fermion mass ma- 
trix, which gives 

trM2/2 = 2e^ (GyG^.^. - 1) - 2g^D''D'' - ^e'^ {g^D^D^Y + Sg^D^D^, . (2.97) 
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The only spin 3/2 particle is the gravitino, therefore one obtains 

trM^/2 = 4M|.e^ = 4m| . (2.98) 

The supertrace formula in supergravity can be written by adding Eqs. (12.951) . (I2.96p . 
dMZD and dZHHD: 

StrM^ = 2{n^- 1) m| - (n^ - 1) g'^D'^D'' + 2/D"Df.* . (2.99) 

Interestingly, when SUSY breaking occurs via F-term condensation (D" = 0), the 
above formula is totally dependent on the gravitino mass. Accordingly, for unbroken 
exact supersymmetry (mg — )■ 0) it gives zero, as is it generally required. The fact 
that in this case the scalar particles are heavier than their fermionic superpartners 
is also relevant for phenomenology. 



2.2 Gravitinos 

According to the formalism which was developed long ago by Dirac and other greats, 
charged fermions correspond to Dirac spinors. The Dirac Lagrangian and the Dirac 
equation constitute the basis for the understanding of the Standard Model of particle 
physics. 

With the advent of supersymmetric theories, which are progressively setting a 
new standard for particle physics, one has to consider also the Majorana fermions. 
They appear as supersymmetric partners of the SM bosons (either scalars or vectors), 
and as such they are represented by neutral (charge auto-conjugated) spinors. In 
particular, the gravitino is the SUSY partner of the spin 2 graviton. Accordingly, 
we begin this section by recalling some basic properties, which are used to introduce 
the gravitino field and perform the calculations in Chapter [3l We remand to the 
Appendix for the notations which are used in this thesis. 

2.2.1 Majorana fermions 

Let us first introduce the charge conjugation operator C for a scalar field (f), namely 

0^(x) = C(j){x) = 7]^(f)*{x) , (2.100) 

where the phase factor is unimodular (Ir/d = 1) and takes the values rj^ = ±1. In 
the chiral representation, C is represented in matrix form as follows, 

C = -^lV ), (2.101) 
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with the Pauh matrix o"^ and the antisymmetric tensor in two components 

e"^ = -6„, = = ( ° J ) • (2.102) 

C interchanges particles and antiparticles. The action on a spinor ip holds as 

^'= = C^^. (2.103) 

ip'^ is called the charge conjugated spinor of ip. In four- component notation, this 
corresponds to 

^ = f it)^^''=( f V (2-104) 



r J " V e 



r=( ;? ) • (2.105) 



One then obtains 

A four- component Majorana spinor ipM satisfies the property [5T] 

= = C^^ , (2.106) 
therefore 1] = C,- In this specific case, we obtain the following: 



The Majorana flip identities are very useful in the calculations of scattering am- 
plitudes. We include them below. If ip and A are two anticommuting Majorana 
spinors, they obey the following 



ipX = Xij, (2.108) 

^75A = A75V', (2.109) 

V-TmA = , (2.110) 

^7/^75^ = A7^75^ , (2.111) 

-ipcrtiiyX = -Xa^yip . (2.112) 

(2.113) 

The above equations imply the useful formula 

V'7^P,.A = -A7^P^^. (2.114) 



2.2. GRAVITINOS 



45 



2.2.2 The free massive gravitino field 

The Lagrangian for the free gravitino can be read off the general SUGRA 

Lagrangian fl2.66p . We remark that it is written in a generic curved space, thus the 
Lorentz indeces (/x, i/, . . . ) in Eg. (12. 66 p are not flat. In the following, we will instead 
work in Minkowski space, consistently with quantization, and use flat Lorentz in- 
dices. This is described in the Appendix, and corresponds to a O*'^ order locally 
supersymmetric system. 

From Eg. (12. 661) . by using the formula (12.931) for the gravitino mass, one finds 

= -^e^''">M757.5p^a - Irn^i^, , (2.115) 

modulo a total divergence. By using the Majorana condition for the gravitino, 
= CipJ^, the above can be rewritten as 

= ^e^'^'^JC'SsT.^p^. + [7^ 7l • (2.116) 

By variation with respect to iptiy this provides with the following field eguation [52] 

e^^'"^757.5p^a + Irn^ b", Y] V'. = , (2.117) 
that is called the Rarita-Schwinger eguation. Operating on Eg. (12. 11 7p yields 

^7^V^M - 7''^V'p = . (2.118) 
Operating 7a 7^, we obtain one more eguation 

2t {dxY^^ - $^x) + (7a7^V'm + 2^a) = . (2.119) 
Multiplication of (I2.119P by 7'^ provides with a third eguation, 

t [Hi^i. - l^H^ + 3mg7^^^ = , (2.120) 

which, by taking into account (12.1180 . gives the first constraint on the massive 
(mg 7^ 0) gravitino field: 

l^^y. = ^- (2.121) 

The second condition follows from substitution of (I2.12ip into (I2.118p . 

3^^^ = ^. (2.122) 
'''This subsection refers mainly to Moroi's PhD thesis [3J and to the references quoted there. 
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By using fl2.12ip in (12.1191) . one finally obtains the Dirac equation for each vector 
component of the gravitino field: 

{0 - mo) i)^ = . (2.123) 

Equations ( I2.12ip . ( 12.122^ and (12.123^ can be indeed solved by constructing a par- 
ticular spinor, that is composite of a wave function u for a spin 1/2 Dirac field, and 
of a polarization vector for a spin 1 field. Working in momentum space, one uses 
the normalization condition 

^(p, s)m(p, s') = 2mQ6ss' , (2.124) 

and the polarization vectors 

e^(p, 1) = ^^(0, cos 6' cos — i sin0, cos 6' sin + z cos 0, — sin 9) , (2.125) 
v2 

e^(p,0) = — (|p|,-Esin^cos0, -Esin^sinc/), -Ecos^), (2.126) 



e^(p, —1) = ^^(0, cos 6* cos (j) + i sine/), cos 6* sin (j) — i cos0, — sin( 
V2 



(2.127) 



for the momentum vector = {E, |p| sin 6 cos 0, |p| sin 6 sin (j), |p| cos 6) of a massive 
particle. Namely, = p^p'^ = m~. These are normalized as 

e;{p,r)e>^{p,r') = -6rr'. (2.128) 

The polarization vectors (I2.125P - (I2.127P obey the constraint 

p%(p,r)=p%(p,r) = 0. (2.129) 

It is straightforward to prove that the wave function for the gravitino, as it has been 
constructed above, satisfies the following equations: 

YMp^^) = 0, (2.130) 
P^Mp^^) = 0, (2.131) 
(^-mg)7/;^(p,A) = 0. (2.132) 

Here '?/'^(p. A) is the gravitino wave function in momentum space. The above con- 
straints imply that in coordinate space, ipf^ ~ e"*^^?/'^ satisfies Eqs. (l2.12T]) . (I2.122p 
and f l2323|) . The normalization of ^/,(p. A) follows from Eqs.f l2.124p and (12:T28|) : 

^,{p,a)r{p,(^') = -2mg5..,. (2.133) 
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The formula for the polarization tensor of a gravitino with mass and momentum 
Pfi [53] can be accordingly written as 



Pf^Pu \ i / Pf^Pe \ / PuP^ \ g £ 

Vt^u - -V ] -nlVt^e [Vun 2^ 7 7^ 



(2.134) 



where the sum is performed over the gravitino helicities ±1/2 and ±3/2. n^j,(p) 
satisfies the following constraints: 

YIIM = ^MY = 0, (2.135) 

p'^n^,(p) = n^,(p)p^ = o, (2.136) 

g) U,,ip) = U^^ip) - mg) = . (2.137) 



m 



In the high energy limit, i.e. when the centre of mass energy is much larger than the 
gravitino mass mg, the spin sum is greatly simplified and acquires an interesting 
structure [53] : 

Il^u{p)^-pV,u + lp^. (2.138) 

G 

One can clearly see here the sum over the ±3 / 2 helicities of the gravitino in the first 
term, while the ±1/2 helicities of the goldstino (the longitudinal degrees of freedom 
of the gravitino) are represented by the second term. Such a splitted form is very 
important to understand the results of Chapter [31 

The solution of Eqs. fl2.135p . f l2.136p . fl2.137p can be expanded in plane waves 
as follows, 

^'^^'^^ = / j2^Wo ^ {'^''''^'^^P' ^^^P^^^^ + e"^^>?(P' >^Hx(^)} ^ (2-139) 

A 

where we have taken into account the fact that ip^ is a Majorana spinor. This can 



be proven also for ijj^, namely ipfj, = ip'^. Clearly, po = y |pp ± In the above 
equation, the time dependence of the coefficient apx{t) is fixed by the Dirac equation 

mm . 



2.2.3 Interactions of the gravitino 



The quantization of the gravitino field is a rather standard procedure, which is dis- 
cussed in the literature [36j, ^37j. Here we do not consider the subject, rather we go 
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straight to the interactions of the gravitino, and write down the Feynman rules which 
are used in Chapter [3l The starting point is the general SUGRA Lagrangian, from 
which we can extract the interactions of the gravitino with the matter multiplets. 
It is clear from Eg. (12. 66 p that this looks like a very complicated task. 

However, the energy scales of the processes which will be considered in the 
following allow to neglect most of the terms in the Lagrangian (12 .66 p . First of all, at 
centre of mass energies ^/s which are much lower than the Planck scale Mp, some 
operators are suppressed by at least a factor ^jMp. Secondly, the gravitinos in 
the scattering processes that are studied in this thesis appear only as external lines 
(on-shell particles). Thus, by invoking the constraint f l2.12ip . we can safely ignore 
the interaction terms which contain o^' V'/i7^- Thus the Lagrangian we use to 
derive the Feynman rules is simply the following [M| [37] : 

(2.140) 

The Feynman rules which follow are reported in Figures flA.7p . (lA.Sp . flA.Op in the 
Appendix. 



2.2.4 Cosmology: the gravitino problem 

So far we have considered the theory of supergravity and how the gravitino emerges 
in this framework. We conclude this chapter by discussing the role of this particle in 
cosmology, which constitutes the motivation of this thesis. The gravitino is relevant 
for several aspects. Not only it has a primary role on theoretical grounds, since 
it is the gauge field of local supersymmetry. Either stable or unstable gravitinos, 
when embedded in cosmological backgrounds, can generate enormous problems to 
theories which are in agreement with observations. This is known in the literature 
as the gravitino problem. 

The production of gravitinos in the early universe occurs through a variety of 
ways, via thermal or non-thermal mechanisms. In the first case, in the reheating 
era at the end of inflation (see Chapter [T]) a thermal bath is generated via the decay 
of the infiaton field. The particles in the plasma scatter off each other, and this 
process generates gravitinos as the final states of 2 — )■ 2 hard scatterings [lU [531 
[55l [56t [571 [58t [59] . It has been shown that an important amount of gravitinos was 
produced right through these reactions [SOI EH]- This is the subject of paper [T], 
where we have found an unexpected effect when the gauge bosons which scatter in 
the primordial bath are massive. 

Non-thermal production is perhaps a more complicated topic. It can follow 
from several sources: from the decay of the next to lightest SUSY particle (NLSP) 
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[en |62l ESI ED , of the moduli fields [651 [Ml EH |68] , of the infiaton [691 [TQl [HI 112] and 
of the SUSY breaking field [73l [Ml [75l [THl [77] . The latter mechanism is discussed in 
the article [2], where it is shown that the reheating temperature can be related by 
the parameters of the model considered, if a right amount of gravitino Dark Matter 
is required. 

The gravitino mass can range in principle from the eV scale up to the 
TeV scale and beyond [H]. It is strongly model dependent, since it is proportional 
to the SUSY breaking scale and therefore to the condensation value of the F field, 
as we have discussed in this chapter. 

In general, gauge mediation predicts the gravitino to be the lightest supersym- 
metric particle, or LSP [16]. Then, if R-parity is conserved, the gravitino is stable 
and it can be a very attractive candidate for Dark Matter [5l]. It can however 
generate problems: since too many gravitinos might overdose the universe, in the 
standard Big Bang cosmology it is set an upper limit of IkeV [55]. This 

constraint can be anyway relaxed if we assume infiation, since it dilutes the initial 
abundance of gravitinoqj- However, the problem persists also in this case. The ther- 
mal scatterings reproduce the particle after reheating. The number density of the 
secondary gravitinos is proportional to the reheating temperature, thus Tr should 
be constrained in order to avoid particle overproduction. 

In other classes of theories, for instance in gravity mediation, the gravitino 
is unstable and if mg is smaller than 10 TeV, it has a lifetime which is usually 
longer than 1 sec. This means that it decays after the beginning of the Big Bang Nu- 
cleosynthesis (BBN) and it generates hadronic and electromagnetic showers which 
can be very energetic P7]. This implies disintegration of the primordial light ele- 
ments [SB] . Since cosmological observations have so far verified the BBN predictions 
to a very high precision, one must impose constraints on the scenario, in order to 
preserve the agreement between theory and observations [791 [HOI [81] . It is therefore 
evident that the gravitino generates several problems in cosmology. In the literature 
lots of effort had been devoted to an extensive study of these issues (see for example, 
[551 [561 [57]). 

In particular, BBN constraints on both unstable and stable gravitinos have 
been recently derived in Ref.[Hn], which contains an updated analysis of the gravitino 
problem in this context. We refer to this paper to complete the above generic 
discussion. For unstable gravitinos, let us consider the decay into a neutralino LSP. 
It is found that the upper bound strongly depends on mg. In order to preserve the 
light elements ^He, ^He, ^Li and D, the upper limit on the reheating temperature 
has a mild dependence on the mass spectrum of the MSSM particles. However, it 



5 Historically, the first bound was calculated to be between 0(1 MeV) and 0(100 GeV) [7S]. 
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strongly depends on the gravitino mass: 

10^ GeV <Tn< 10^° GeV if 300 GeV < < 30 TeV . (2.141) 

The above range for Tr is interesting, since it provides with a very stringent con- 
straint on theories of thermal leptogenesis, which actually require a reheating tem- 
perature that is at least of the order 10^ GeV plEallHlEHllHniEIllHHlIM]. 

In the case of stable LSP gravitinos, the bounds depend on which particle is 
the NLSP (the authors have considered the bino, the stau and the sneutrino). Very 
stringent constraints have been found if the bino or the stau are the NLSP. Namely, 
mg 10 GeV is excluded if the mass of the NLSP is lighter than 1 TeV. If the 
sneutrino is the lightest supersymmetric particle, the BBN constraints are sensibly 
weaker since the sneutrino decays mainly into gravitino and neutrino (i.e. weakly 
interacting particles). In any case, the constraints are generally very restrictive. The 
decay of the NLSP alone cannot provide with the correct amount of gravitino Dark 
Matter. Therefore the largest contribution should be produced by other mechanisms, 
for instance through thermal scatterings or by decay of the scalar condensate. 

Both of these possibilities are investigated in the research papers [1] and [2], 
which are discussed respectively in Chapters [3] and H] of this thesis. 



Chapter 3 



Gauge boson scattering and 
gravitinos 

In this chapter, the phenomenology of the Standard Model and of the MSSM will 
be discussed. Topics like gauge invariance and the MSSM mass spectrum are closely 
related to the results of paper [1], that concerns the scattering of two massive W 
bosons 

W + W^-^W^ + G, + W- -^x^^ + G , (3.1) 

with a gravitino and a gaugino in the final state, both in the gauge and matter 
eigenstates. In the next section we analyse the massless limit of the above process, 
namely the gluon scattering in supergravity. 

3.1 Scattering of gauge bosons in supergravity 
3.1.1 The massless limit: gluon scattering 

As a starting point, we consider in analogy with Ref. [53] the production of gravitinos 
through a process in QCD. Namely, we study the scattering of two gluons g*^ and 
g^p with a gravitino G and a gaugino (gluino) g'^ in the final state. This is process 
[A] in [53], 

g^ + g'-^r + G, (3.2) 

where a, 6, c are indices of the SU{?>)c algebra. 

The reaction f l3.2p corresponds to the four Feynman diagrams in Fig. (l??p . 
Three of them contain the channels s, t, m, while the forth is a contact diagram, with 
a point-like supergravity interaction [M] . 
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The matrix element Ai of the process f l3.2p can be written as the sum of the four 
subamphtudes 

M = Ms + Mt + Mu + , (3.3) 



where s, t and u are the Mandelstam variables |90]. By using the Feynman rules in 
Appendix IA.2t it is easy to write the four amplitudes as follows: 

Ms = [v'^^ik - k'Y + 7]^'' {2k' + kT - t''{2k + k'f] X 



+ r,7.] l'v';,{p')] el{k)e%k') , (3.4) 



= -f^{^;(p)[r,7"]7^(^-/ + m,)7X^(pO} '°^^!t^r ^ (3-6) 
= -|^{^^(p)[7",7^]7^<'(p')}e^(A;)4(^')- (3.7) 

In the above equations, Mp = (SuGn)'^^'^ = 2.43 x 10^^ GeV is the reduced Planck 
mass, and gs and fate are, respectively, the strong coupling constant and the struc- 
ture constants of the group SU(3)c. The spinors are ip^^ip), namely the wave function 
of a gravitino with four-momentum p and helicity state s, which is represented by a 
four-component Majorana spinor, with spin 3/2 and Lorentz index /x. v'^,{p') is the 
wave function of an antiwino with four- momentum p', helicity state s' and SU(2)l 
index c in the final state. It is represented by a spin 1/2 Majorana spinor (see 
Section EXID. 

In the following, ipl^{p) and v'^,{p') are replaced respectively by ipf^ and v"^ for 
a clearer notation. The wave functions of the gluons, e^(/c) and e^^^k'), are the 
polarization vectors of two bosons with, respectively, four-momenta k and k', SU(3)c 
indices a and b and Lorentz indices a and /3. The kinematics of the process defines 
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the Mandelstam variables as follows, 

s = {k + k' f , 

t = {p'- k'f , (3.8) 
u= ip' — kf . 

The conservation of energy- momentum is written as [90] : 

s + t + u = = m? + m~ . (3.9) 

i 

Ward identities and connection to unitarity 

Before calculating the squared amplitude relative to the process (13. 2p . let us recall 
the generalized Ward identities for Yang-Mills theories [901 E] • In the non-abelian 
case, they are called Slavnov- Taylor identities [92]. We basically follow the book by 
Cheng and Lee [91], where a detailed calculation can be found. This section will 
provide only with the comments which are relevant to our results. 

The Ward identities relate different Green's functions, namely these reflect the 
symmetries of the theory such as gauge invariance. For instance, in electromagnetism 
these imply one single condition, namely that particles with the same bare (electric) 
charges have the same renormalized charges. In non-abelian gauge theories the 
number of these conditions is much larger. Moreover, the Ward identities guarantee 
that all the unphysical singularities eventually disappear in the physical amplitude. 

Let us consider an SU(2) theory with fermions / in a doublet representation. 
The requirement that the scattering matrix S must be unitary, 

SS^ = S^S=1 ^ J2^-i^bi = Sab, (3.10) 

I 

by virtue of the relation 

Sab = Sab + i{2n f5'^{pa - Pb)Tab , (3.11) 

can be written in terms of the scattering amplitude Tab as follows: 

ImT„, = l5^T,,T,*(27r)^54(p„ -p,). (3.12) 

The unitarity of the 5'-matrix thus implies certain conditions which any scattering 
amplitude has to satisfy. In particular, the imaginary part of T"^ connects the initial 
and final states as a consequence of energy conservation. 
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Figure 3.2: The Cutkosky rule applied to fermion-antifermion scattering. 



Suppose now to study a fermion-antifermion scattering // — )■ // in this basic 
SU{2) theory, with off-shell gauge bosons in the intermediate states [93]. In order 
to recall the property (I3.12p . one can use the Cutkosky rule [Q^ to calculate the 
imaginary part of the scattering amplitude. One then replaces the boson propagators 
by their imaginary parts, which are multiplied by the on-shell scattering amplitudes 
T{ff — )■ AA) and T*[AA — > //). This is equivalent to making a cut in the original 
diagram, as in Fig. (3.2). 



It can be shown [91] that at the lowest order the unitarity condition for this scattering 
reduces to 

j dp KT^t*^''"^"'' - 2S^'S'''*) = j dpT^lT^l*P^''{ki)P''^{k2) , (3.13) 

where P'^°'{ki) and P'^^{k2) are the spin sums of the gauge bosons: 

p/^"(A;i) = ^e(")'^(A:i)e('')"*(A:i), (3.14) 

r 

p-^(k2) = J2e^'->{k2)e^'-^^*{k2). (3.15) 

r 

In the unitarity condition Eq. fl3.13p . dp is a volume element in the phase space of 
the two bosons. T^^ corresponds to // — >■ A^j^A^ and 5'"'' is the amplitude of the 
process // — )■ c"'c^\ where c are the ghost fields. In the left hand side of (13.131) . 
these give an explicit contribution to the scattering amplitude. In the right hand 
side, they appear instead as the longitudinal degrees of freedom of the gauge bosons 
in the above spin sums. 

In other words, the spin sum of the massless gauge boson in non abelian 
theories cannot just be r/^j, (as in QED), but it must contain also the unphysical 
ghosts in the form of longitudinal polarizations of the gauge bosons. 
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This is our starting point. In the case of gluon scattering in supergravity, 
we encounter exactly the same form of equation (13.131) and verify the unitarity 
condition. In contrast, in Section 13.1.21 we show that in the case of massive W 
bosons, the longitudinal degrees of freedom in the r.h.s. of Eg. (13. 131) are missing, 
by virtue of the SUGRA vertex in the s channel. 

Squared amplitude of the gluon scattering 

Let us now return to the gluon scattering in supergravity, by keeping in mind the 
discussion in the previous section. The sum of the four contributions Eqs. (13.41) . 
( 13. 5p . (13. 6 p and ( 13. 7p is the matrix element A^, that can be written as follows: 



The Ward identities for the above amplitude assume the following form |91] . 



where ka and fc^ are the four momenta of the two on-shell gauge bosons. Sab and Tab 
are the Lorentz scalars which have been defined in the previous section. According 
to the above discussion, here we call them ghost amphtudes. By using equations 
(13.171) . we can write the squared amplitude, summed over the polarizations of the 
initial and final states, in a very meaningful way. A straightforward calculation 
gives: 




(3.16) 




(3.17) 



spin spin 



(3.18) 



spin spin spin 



By using Eq l3.16p . the scalar amplitude can be rewritten as 




(3.19) 
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therefore fl3.18p can be recast in a very simple form: 

J2\M\' = J2\K,r-2j2\Sa,\'- (3.20) 



spin spin spin 



The above formula follows from the polarization sum of a massless vector boson 
with momentum k and helicity state r [91], 

E ^t\k)ei^'^*(k) = - (^va. - ^"^- + ^-^" ^ 6^^^T + L, (3.21) 

where the two different contributions of the transverse and longitudinal components 
are labeled by T and L. Clearly, J2spin l^abl"^ comes from T = —rjau, while L gives 
~'^Ylspin\^ab\'^- '^^^ polarizatiou sums of the fermions result as follows. For an 
anti-gaugino (anti-gluino) [42j with momentum p', group index c, helicity state s' 



and mass nig, we have 



5^<,(p')^r,(p') = (/-m,)5'=". (3.22) 



For a gravitino with momentum p, helicity state s and mass mg, the general expres- 
sion for the full projector fl2.134p has been given in Section 12.2.21 as 



+ mg) 



P^LPu \ 1 / P^^Pe \ / PuP^ \ e £ 

"^G H K \ K 



(3.23) 
(3.24) 



One can simplify the above object by assuming that the centre of mass energy is 
much larger than the gravitino mass. In this case, the full gravitino projector fl2.134p 
reduces to 153): 



U,,{p)^-fV,. + lf^- (3.25) 

G 

This approximation is consistent with a large class of cosmological scenarios with a 
reheating temperature that is well above the electroweak threshold. This corre- 
sponds to Refs.[53] and [59]. Even though in models with low Tr this is not always 
true, for the gluon scattering at hand we make the above assumptior]0. The reason 
for using equation (13.25^ is that the contributions from the gravitino polarizations 
±3/2 and ±1/2 are clearly distinct. 



^It can be shown that the fuU gravitino projector (|3.24l) does not give a quahtatively different 
squared ampHtude. 
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The kinematics is derived by taking into account Eq. fl3.8p . 

2kp - 



e = 


= , 


2pp' 




m|, 


2k'p 




m?, 


2kp' 


2kk' 


= s, 


2k' p 



[s — m? — m~) 



(m? - t) 



(3.26) 



The final expression for the squared matrix element of the scattering 

summed over the spins of the initial and final states, and averaged over the initial 
states, can be accordingly written as: 

.2 



spin spin 



gM 



abc 



Ml, 



1 + 



3m~ 

G / 



s + 2t + 2— 
s 



(3.27) 



This agrees with the result obtained in [53]. In the next section we perform the 
same analysis for on-shell massive bosons and see that things change non trivially, 
despite what is generally expected. 



3.1.2 Scattering of massive gauge bosons 

Here we calculate the cross section of the WW scattering into gravitino and gaugino, 
and compare it to the results of the previous section. Since the supergravity vertices 
are universal, one expects that the scattering amplitude of 

Wik) + W\k') — > W%p') + G{p) , (3.28) 

where a, b, c are now indices of the SU(2)l algebra and is a wino, should not 
be formally different from Eq.f l3.16|) . The Feynman diagrams of the above process, 
whose structure reflects this universality, are drawn in Fig. (l3.3p . The four corre- 
sponding amplitudes are practically the same as in Eqs.f l3.4p . (13. 5p . (13. 6p . (13. 7p : 

M, = ^ [r]''^{k - k'y + r]^''{2k' + k)" - r]'"'{2k + k'f] x 



s — m 



w 



Mt = V{^^(P)[^'7l7nr-/ + m^)7V'(p')}e^(^)4(^'), 

^ w 

Mx = -f^{^^(p)[7",7^]7^<'(p)}e:(^)4(^')- (3.29) 
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G w 




Figure 3.3: The four diagrams which contribute to W + — > + G. 



These can be rearranged in a single matrix element, which takes the form 



s niy^ t m~ 



u — m 



w 



(3.30) 



The explicit dependence of the polarization vector of the W on the momentum has 
been omitted. As it was done before, the above amplitude is recast by using the 
Ward identities Eq.f l3Tfl) . 



k'pM:,^ = Tabk^ 



(3.31) 



In this case. Sab and Tab are no longer due to the ghosts, but to the physical degrees 
of freedom which are carried by the Goldstone bosons. Accordingly, we will regard 
them as Goldstone amplitudes. These have the following expression. 



^ _ geabc I 1^'] 7^^' . _ ^ 



4Mf 



s — m 



(3.32) 



w 



As expected, the above differs from fl3.19p only trivially, by the W boson mass niw 
in the denominator. The spin-summed square of the matrix element in Eq. fl3.3Qp 
can be now written as: 



(3.33) 



spin 



spin 



spin 
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By using the full gravitino projector f l2.134p and averaging over the initial states, 
we get 



\M\^ = 4 



g2 1 ^abc 1 2 



5^'^"' " 9M2 

spin ^ 





(3.34) 

rrii can be either m^/ or m^, and / is a function of dimension [mass]'^. In the 
previous section we have obtained for SUSY QCD the following, 

E - i E l-^r - ^ (l + ^) + 2« + 2^) , (3.35) 

spin spin \ G / ' 



Therefore the relevant difference with respect to Eg. (13.341) is found in the term 

-t{s + 1) 



Cjr 



(3.36) 



that contributes to the differential and to the total cross section of the process as 
follows: 

\dt),^ UtiMI 3m| ^ M>| • ^^-^^^ 

Clearly, the above result implies violation of unitarity of the scattering matrix for 
sufficiently high centre of mass energies, since the cross section grows with energy. 
Fortunately, the structure of (I3.36[i n makes the unitarity breaking energy coincide 
precisely with the supersymmetry breaking scale (see Chapter [2]). Since supergravity 
is the low energy limit of some GUT theory yet to be discovered, in principle the 
above result does not affect observations. Unitarity is guaranteed up to the scale 
of validity of this effective theory. For instance, since the SUSY breaking energy in 
supergravity corresponds to |35j 



/^susY = ^{F) = iV^rnQMp)'/^ , (3.38) 

we find that unitarity is broken e.g. at y/s ^ 0(10^'' GeV) if ^ 0{1 MeV). 

Let us now summarise our results. After calculating the gluon scattering in 
supergravity (13. 2p . which provided an expected result, we have considered its exten- 
sion to massive gauge bosons. We have found that the cross section for the WW 
scattering contains an extra term which breaks unitarity at the SUSY breaking scale. 
Now we would like to understand this unexpected result. 

To begin with, we observe that Eq. (l3.37|) has some similarity with the case of 
WW scattering in the Standard Model [90j. The Goldstone equivalence theorem, 



^Namely, the fact that the Planek mass appears in the denominator. 
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that will be discussed in the next section, states that the longitudinal polarizations 
of the massive on-shell W are crucial. We expect that these will be relevant for the 
case at hand as well. 

Consider the squared amplitude written as f l3.33p . The main difference with 
respect to the gluon scattering (13.21) . 

spin spin spin spin spin spin 

(3.39) 

is found in the different contribution of the scalar amplitudes. It can be shown that 
M"ff is formally the sam (i as Eq.f l330|) and that the masses of the W bosons in the 
kinematics do not give any contribution which would cancel fl3.36p . The sum over 
the polarizations of a massless gauge boson with four- momentum k is Eq. (l3.2ip : 

E C^'\k)el^'^^^*ik) = - {^a. - ^"^- + ^-^" ^ 6'^'^T + L. (3.40) 

On the other hand, one obtains Eq. fl3.33p from the projector of a W boson with 
momentum k, mass mw and helicity r: 

J2e^Hk)e^'>{k) = - - ^) 5"' ^ T + L. (3.41) 

We have written (I3.33P and fl3.39p in such a way that the two distinct contributions of 
the above projectors stay separated. Thus l^^ff P comes from the transverse 

polarizations riay (i.e. from T). The longitudinal degrees of freedom^ give instead 
Sspm I'S'abP- As a cross check, by substituting the spin sums (I3.4ip directly into 
the square of fl3.30p . we obtain the same result. In Table fl??p we compare the 
divergences generated by the distinct terms T and L, for both massless and massive 
W. These lead respectively to Eqs. (13.270 and fl3.34p . 

It is now evident that fl3.36p is somehow related to the longitudinal degrees of 
freedom of the W, namely to its mass. We try to find a theoretical explanation for 
Eq. (l3.34p in the next subsection. 

On the mass of gauge bosons and unitarity violation at high energies 

We have seen so far that the gluon scattering into gravitinos provides with an ex- 
pected result fl3.27p . in which the squared amplitude factorizes a mass splitting that 

■^Modulo constant factors such as the couphngs and the structure constants. 
^It is anyway important to point out that the gluons do not have any proper longitudinal 
polarizations, since these are massless. Such degrees of freedom are the Faddeev-Popov ghosts. 
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Figure 3.4: The contributions of the spin sums ( I3.2ip and ( I3.4ip to Jospin l-^P- 



vanishes in the limit of exact supersymmetry. On the other hand, if the gluons are 
replaced by W bosons in the broken phase, a new term that breaks unitarity at a 
certain energy scale appears. Eg. (13. 36 p . 

Since the interactions of the gravitino with the gluon and the W are exactly 
the same (modulo constant factors), we concentrate on the longitudinal degrees of 
freedom of the W. These are carried by the Goldstone bosons after the spontaneous 
breaking of the gauge symmetry and are connected to the non vanishing mass of 
the W. In the following, we show that their behaviour is non trivial at both low 
and high energies. In other words, the divergence (I3.36P appears not only in the 
electroweak phase, but also at high energies. 

We consider the amplitude of the s-channel, with an off-shell W. It contains 
the product: 



s — m 



w 



Vau + (e - 1) 



{k + k')^{k + k') 



{r,{p)[i(^ + f(^',r]YvUp')}- (3.42) 



s — ^mf 

The W boson propagator in the gauge can be written as follows: 
< A^{q)A^{-q) >=< A^{k + k')A^{-k - k') >= 



s — m 



w 



Vau + (e - 1 



{k + k')^{k + k') 



s — 



2 

w 



(3.43) 



This is a completely general formula, that reduces to the usual propagators in the 
Landau or Feynman gauge by choosing a particular value (or limit) for C,- In partic- 
ular, the unitarity gauge corresponds to ^ -^00: 



<A^{k + k')A^{-k-k') >-- 



m 



w 



{k + k'),{k + k'\ 



m 



w 



(3.44) 
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Here we see that, for both massless and massive gauge bosons, the part of the 
propagator which takes into account the longitudinal degrees of freedom of the 
particle. 



-I 



s — m 



w 



(3.45) 



is cancelled identically when contracted with the SUGRA vertex 



V{G, Al, W^) = --^ + r, . (3.46) 

This happens by virtue of the commutator. The result is thus the same in the 
Feynman (^ = 1), Landau (^ = 0) and unitarity gauge (^ — )■ oo). It is clear that 
this is a general result in supergravity. 

It is well known from the Standard Model that for massless gauge bosons, the 
longitudinal degrees of freedom are unphysical and cancel during the calculation. In 
contrast, if the gauge bosons are massive such cancellations should not occur, as the 
longitudinal polarizations of the W carry physical degrees of freedom. This might 
be responsible for what happens when in the squared amplitude we find terms such 
as 



t{s + t) s 



CTtot ~ — 2" • (3-47) 



— ?,Mlm% Mlm% 

spin ^ G ^ G 

In fact, the Feynman and the unitarity gauge are not equivalent. If the propagator 
is expressed in the Feynman gauge as [90] 



< AJk + k')AJ-k - k') >= -i , (3.48) 

s - 

one must recover the missing longitudinal modes by adding further diagrams with 
the exchange of a scalar. In the Standard Model, once the gauge symmetry is 
broken, the unphysical Goldstone boson is eaten by a massive gauge boson. The 
Goldstone controls the amplitude for emission or absorption of the gauge boson in 
its longitudinal polarization state |90] . 

Consider now the degrees of freedom which are involved in the Higgs mecha- 
nism. We see that a massless gauge boson (two transverse polarizations) combines 
with a scalar Goldstone (one polarization) to produce a massive vector boson, with 
three polarizations. When the boson is at rest, the three polarization states are 
completely equivalent. At relativistic velocities the longitudinal degree of freedom 
"decouples" instead as a single Goldstone boson, see Fig. (3.5). This effect is called 
the Goldstone boson equivalence theorem and was first proven by Cornwall, Levin, 
Tiktopoulos and Vayonakis [9S]. Let us therefore consider a massive W boson. When 
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Figure 3.5: Visualization of the Goldstone boson equivalence theorem. At high 
energy, the amplitude with absorption or emission of a massive gauge boson is equal 
to the amplitude with absorption or emission of the Goldstone scalar which was 
eaten by the gauge boson. 



it is at rest, the momentum is k'^ = {mw, 0, 0, 0) and the polarization vector is a 
linear combination of three orthogonal unit vectors: 

(0,1,0,0), (0,0,1,0), (0,0,0,1). (3.49) 

The first two give the transverse polarizations. Boosting the boson along the z-axis 
changes the momentum to /c'^ = (ii^k, 0, 0, k). The polarization vectors are three unit 
vectors and satisfy 

e'^fc^ = 0, e-e = -l. (3.50) 

Clearly, since the boost is along the third direction, the transverse vectors do not 
change. Instead, the third one corresponds to the longitudinal degree of freedom, 

(^,0,0,^) , (3.51) 



which in the approximation A; — > oo can be recast, component by component, as 

The components of are growing with k, in agreement with the general relations 
eL ■ k = and k ■ k = m^. 

By keeping in mind this last formula, we can now discuss the high energy 
behaviour of the result 03.341) . and prove that the longitudinal polarizations of the 
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two on-shell W bosons in the scattering 



W + G. 



(3.53) 



become strongly interacting at high centre of mass energies. This imphes that the 
term 

ti-s + t) 



3Mlml 



2 ' 



(3.54) 



G 



holds also above the electroweak threshold, namely at any energy scale. 

In supersymmetric theories with unbroken gauge symmetry, the supercurrent 
is conserved in the SUSY limit. Namely, the right hand side of 

p^S, = A(M) , (3.55) 

where is the four momentum of the gravitino and A(M) is the mass splitting 
of the supermultiplet, vanishes when all the masses are set to zero. However, in 
the electroweak theory the W mass cannot be set to vanish. In other words, we 
expect that in the process fl3.ip the longitudinal modes of the W become strongly 
interacting in the limit mw — ^ 0. Accordingly, in this case the statement f l3.55p 
does no longer hold. We prove this claim immediately. The matrix element can be 
written as 

m = Mp)s^, 

therefore the squared amplitude can be recast as 

\m\' = s''Mp)-Mp)S'- 

When summed over the gravitino polarizations, this becomes 



(3.56) 



ip—spin 



i^v{p) ■ i^^,{p) 

jp—spin 



(3.57) 



Since a term with the gravitino mass squared in the denominator can appear only 
from the longitudinal spin 1/2 components of G in either f l3.24p or f l3.25p . one can 
replace the gravitino projector with the second term of Eq. fl3.25p . that is: 



■0— spin 



OC 



2 

m~ 
G , 



S" 



m- 



G 



(3.58) 



Therefore, it is always proportional to Pf_iS^ = — )■ p^), that is to the matrix 
element with the replacement ip^ ^ Pp.- 
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It is now easy to prove that the scalar density Pfj,S'^ is not proportional to any 
SUSY breaking parameter, and it does not vanish in the SUSY limit. By explicit 
calculation, it can be shown that Eqs. (13.291) provide with the following: 



5af + 0{m) 



4Mp (s - m2^)(t - m|;)(-s - t + Iw?^ + ml) 



e^(A:)e^(A:'). (3.59) 



The tensor density has a rather nontrivial expression: 



= 2m^ [s2(-7°7^|^' - 2p"'7°) + t^{ri'^^^ - r]"^^^' - 2A;^7" + 2A;'"7^) + 
+st(-r/"^|^ - f]""^^' - 2A:^7" + 2k"'-f'^ + -f"-f^^ - 7°7^|^' - 2p'"7^ - 2p'^7")] . 

(3.60) 

At this point, one would be tempted to neglect the terms proportional to Tfiy^ as 
mw — 0, and the claim of the conservation of the supercurrent would immediately 
follow. However, if we now recall Eg. (13. 52 p and substitute it into (13.591) . the result 
is (modulo constant factors): 



rriQ Mpin^ \"^g 

where now m is really any mass. We have thus seen that in the scattering of 
two massive gauge bosons in SUGRA, the supercurrent is not conserved in the 
SUSY limit. This happens since the non vanishing term in the right hand side of 
(13.611) violates the conservation of the supercurrent. It is also easy to show that the 
corrections O (m/E-k) provide with terms proportional to m^, which are absorbed 
into O (m^/mg). Therefore the proof given above is valid also at low energies. 

To summarise, in this section it has been shown that the structure of the 
supergravity vertex cancels the longitudinal degrees of freedom of any propagating 
gauge boson. If the boson is massive, terms with bad high energy behaviour are 
generated. With an argument based on the supercurrent, we have proven that the 
above result holds both at low and high scales. At high centre of mass scales, the 
longitudinal degrees of freedom of the W become strongly interacting and violate 
the unitarity of the theory. 



^We have used equations like k"ea{k) — k'^ep{k') = 0, f'v'^{p') = m^v'^{p') = 0{m), and 



m 



w 
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Figure 3.6: The annihilation diagrams with a propagating Higgs. 



3.2 Gauge eigenstates and the Higgs boson 



In the previous sections we have understood that the terms which lead to violation 
of unitarity in the cross section of the process 



W' + G. 



(3.62) 



are related to the missing longitudinal polarizations of the off-shell W in the s- 
channel. Accordingly, if we want to solve the problem by cancelling the divergences, 
we must recover those missing degrees of freedom in some way. In the case of 
( 13. ip . the only possibility is an annihilation channel with an off-shell scalar. We can 
anyway see immediately that in the basis of gauge eigenstates there might be some 
formal difficulties. 

First of all, in order to be consistent with the process at hand, the scalar 0^ 
has to be included in the adjoint representation of the gauge group. However, in 
Supergravity the scalars are put into the fundamental representation. In this case 
anyway, the gauge symmetry is broken, so the particle multiplets can pass from a 
representation to the other, by virtue of the mixing of gauge eigenstates. We consider 
the dominant couplings of the gravitino with the supercurrent in the supergravity 
Lagrangian fl2.66p . 



V2Mp 



8M, 



up 5 



(3.63) 

where Xl r~ Pl,rX' = 1- 1^) ■ The coupling of the gravitino to a scalar and 
a gaugino does not exist in the above equation. We can anyway use the method 
of mass insertions, and invoke the mixing of the gaugino with a chiral fermion that 
couples to the gravitino and to a scalar. 

In contrast, the supergravity couplings in f l3.29p regard Majorana fermions 
which do not have definite helicity. Namely, in the interference of the Higgs scalars 
with the original amplitude f l3.3Up that is now labelled with M^^^ , some polarizations 
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which might in principle contribute to the process are cancelled. Accordingly, we 
expect that the additional diagrams in figure (13. 6p would not cancel the divergences 
(I3.36p . In the following, we consider the particle content of the MSSM [12], that is 
summarised in the table at the end of this section. It contains two neutral Higgs 
bosons and i?^, which are suitable to our proces^. These break the gauge 
symmetry by taking the vacuum expectation values Vi and V2, 

where Vi ^ V2. We label with M^^^ the amplitude which results from the two 
distinct diagrams with the exchange of a Higgs, as in Fig. (l3.6p . 

The vertex of the Higgs with the W bosons is non trivial either, since we 
consider a doublet of scalar s instead of a singlet. Accordingly, this is not the same 
coupling of the Standard Model. This leads to the following expression for the vertex 

V{HtW:wl) = '^C^^rTf'e'^^^ C], = v.Z^ + v^Zl . (3.65) 

Zji is the rotation matrix transforming the gauge eigenstates iff and if| into the 
mass eigenstates H and /i, namely the heavy and light Higgs in the MSSM: 

[ f Uf -sina\ / m\ 

\ h ) \^ sma cos a ) \ ti^ ) 

We are actually interested only in the gauge eigenstates, therefore we impose a = 
in (I3.65P to obtain the following: 

V{H^W^W^) = ^-^e^^'^viT]''^ , (3.67) 

V{H^iW^Wl) = !|-e"'"=T;2r/"^ . (3.68) 

Recalling the mass eigenstates of the MSSM, we remember that the higgsinos and 
the winos mix with each other. This justifies the mass insertions. Physically, the 
higgsino hi, which interacts with the scalar if? and with the gravitino G, at a certain 
point becomes the wino. The mass insertion keeps trace of this mechanism. The 
mixing factors can be found in the neutralino mass matrix [42] . 



A2 = 

^Throughout this section these are cahed 



+ Y for hi, (3.69) 
forh2, (3.70) 
Hf and HS, where c is an index of SU(2)l. 
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where g is the couphng constant of SU(2)l. Moreover, the higgsino which propagates 
in the external leg of each diagram provides with the factor 

^ ^ (^ = 1,2). (3.71) 



The above approximation will give a clearer form of the amplitude, and it can 
be safely assumed since the higgsino mass does not enter the kinematics and is not 
relevant for the final result. We thus obtain the amplitudes corresponding to the 
neutral Higgs bosons of the MSSM in the basis of gauge eigenstates, 

M(^f) = ( ^1 ^.(1^ + ^)7^1-75^ , ^3^72) 

MiH,) = ( M±lV^L^l^^«.,.(fc),^(fc') . (3.73) 

By recalling the fundamental relation 



^w = ^iv! + vl), (3.74) 



this can be recast into a single amplitude: 



g^abc 



8V2Mp 

" + - mj,M - Amlmi:\ \ ^.(^ + (1 ' 75)^^^ _ 

(3.75) 



4 



Now we can go back to the original calculation. It can be proven that the anomaly 
in fl3.36p is not cancelled by adding the above expression to fl3.30p . The dominant 
divergences result as: 



„2 1 _abc 1 2 



X 



^t{s + t) 1 / _ s2(7s + 3t) _ 64s(s + 2t) 3s 

3m| 64 x 48 ym^m| m|;(s + t) ^m^m^ 




(3.76) 



The interferences of M^^^ with M^^^ give only the term with \/2. This is expected 
from our previous discussion, since the Higgs vertex contains the helicity projectors 
which cancel a number of contributions. However, by iterating the method of mass 
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+ 





Figure 3.7: The t— and u— channel diagrams with the mass insertions. 



insertions, one may add further diagrams, developing a perturbation theory around 
the mass ratio {my^jm^\ At the first order, the matrix element M^^^ is given by 



(3.77) 



where Mt^ , and , refer to the t and u channel. The neutral higgsinos h\ and 

mix both with the wino W ^ as it is shown in figure (13. 7p . Therefore there are two 
amplitudes which can be summed to give 



1 g^ahc ( m'^ 



tl,2 



4^2 Mp \m 



"w 



1 



t — m~ it — m~)t 



for the diagrams with the t-channel. Regarding the u-channels, we get: 



(3.78) 



1 ga-abc ( rn^ 



4a/2 Mp \m 



"w 



m 



w 



if - f) 



u — m~ iu — m~)u 



(3.79) 

The amplitude at second order, M*^^\ contains the double mass insertion with a 
gaugino in the external leg. This corresponds to the diagrams in figure (13.81) . where 
i = 1,2. The corresponding amplitudes are the following, 

^ T^::^ (^M + ^ . 7.] 7^^^) + {^M 7l 7^(r - / + m^)l^v^] 



s — m 



w 



u — m 



^ {^,[r,7^]7^(^ -/ + ^w)Tv^] - ^^l^l'-.l^h'v'Aelel . (3.80) 



w 



The above is nothing but Eg. (13.301) . multiplied by the factor m^rlm—. The total 
amplitude can be finally written as the sum of the four contributions: 



(3.81) 
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v'' Wji 




Figure 3.8: The diagrams with two mass insertions in the external legs. 

The squares and the interference of M^^^ and M^^^ give subdominant terms which 
are proportional to for instance 



t{s + t) 



mi 



w , 



3m% 



(3.82) 



Thus the leading divergences appear only in Eq.f l3.76|) and in 



spin 



48M|. 



(3.83) 



Remarkably, the above factor cancels exactly the corresponding term in Eq.f l3.76|) . 
This means that the contribution of the amplitude M^^^ is suppressed, and that 
the Higgs diagrams contribute only with the square Xlspm then 
conclude that, as expected, the neutral Higgs doublet of the MSSM gauge eigenstates 
do not restore unitarity. 



3.3 WW scattering in the basis of MSSM mass 
eigenstates 



We have seen in the previous sections that the scattering 

Wik) + W^k') W'{p') + G{p) , (3.84) 

in the case of massive W bosons, provides terms like (13.361) in the squared amplitude. 
Such quantities can violate the unitarity of the scattering matrix above a certain 
scale. The basis of gauge eigenstates is useful from a theoretical viewpoint, since it 
separates the contributions of the distinct gauge groups. The calculations turn out 
to be also simpler. However, we have seen that in this framework the anomalous 
terms are not cancelled. 
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Ordinary particles 


Supersymmetric partners 




Gauge eigen- Name 
states 


Mass eigen- Name 
states 


q=u,d,s,c,b,t 
i = e, ii,T 
9 

W± 

Hd 

7 

HO 

u 

B 


qL, qn squark 
(-lAr slepton 
g gluino 

wino 

H+ higgsino 
higgsino 

7 photino 
Z° zino 
H° higgsino 
higgsino 
wino 

B bino 


gi,g2 squark 
£i,£2 slepton 
g gluino 

Xi2 charginos 
X° neutralinos 



Figure 3.9: The particle spectrum of the MSSM. 



Moreover, the observation of particles at colliders is based on the electroweak 
gauge eigenstates, while the mass eigenstates are the propagating states, when the 
observation is interpreted within the theory. In the Standard Model, for an unbroken 
gauge symmetry such as SU(3)c, there is no mixing. On the other hand, when the 
gauge symmetry is spontaneously broken, the mass matrices are non diagonal. The 
eigenstates of such matrices are new physical particles, obtained after a suitable 
rotation that follows a change of basis in the field space. 

This happens naturally also in supersymmetric theories, and in particular in 
the MSSmIiI. In the basis of mass eigenstates, the particle content includes the four 
neutralinos x° {i = 1, ...,4) and the four charginos xf (j = l^S). This is reported 
in the table in Fig. (3.9). The process 

W+ + W- -^x'l + G , (3.85) 

is the counterpart of Eq. (13. ip and corresponds to the set of diagrams in figure (I3.10p . 
X° is a neutralino with i running from 1 to 4. This is clearly a physical process that 
can be observed in experiments such as the LHC [57]. Since the amplitudes for each 
single Feynman diagram are complicated, these will not be included here unless 



detailed discussion of the mass eigenstates of the MSSM can be found in Ref. |42] . 
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Figure 3.10: WW scattering in supergravity in the basis of mass eigenstates. 



necessary. Before making any calculations, one might infer about the possible result 
by looking at the Feynman diagrams in Fig. fl3.10p and at the supergravity vertices 
in Chapter |2l The vertices now factorize the mixing factors of the mass eigenstates, 
but do not change analytically. This means that the cancellation of the longitudinal 
modes in the Z° boson propagator still occurs, thus we expect that also in this case 
the longitudinal polarizations of the on-shell W become strongly interacting in the 
high energy limit. 

As already discussed, the higgsinos and the electroweak gauginos mix with each 
other by effect of the electroweak gau^ symmetry breaking. The neutral higgsinos 
and H^, and the neutral gauginos W'^ and B mix and form four mass eigenstates, 
the neutralinos. The charged higgsinos if^ and H^, and the winos W^ form instead 
the four charginos. 

Let us start with the neutralinos. In two component notation, following 
and [H], we can choose the basis of neutral mass eigenstates 



(3.86) 



where A' and A3 are linear combinations of A^ and A^. One can then write the mass 
matrix as follows: 



Y 



f M' 

M 
— m2 sin 6*^ sin sin 6'^, cos ^vk 

y cos 6*1, sin 6'vi/ — cos 6^^, cos 



—rriz sin 6^ sin 6^ cos 9^ sin 9w \ 
ruz sin 9y cos 9w —mz cos 9^ cos 9w 
-/i 
-/i J 

(3.87) 
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where tanOy = Vi/v2. The mass eigenstates are defined as 

X° = iV,,V;°, ^ = l,---,4, (3.88) 

where is a unitary matrix so that one obtains, for the diagonal matrix Nd with 
nonnegative entries, 

N*YN-^ = Nd . (3.89) 

Since the Feynman diagrams in Figs. (13. Ill) and (13.131) contain propagating photons, 
Z° and the heavy and hght neutral Higgs bosons H and h, we choose the basis 

= {-z\^,-iK,\H,\h), (3.90) 

that contains the photino 7, the zino z and the higgsinos H and h as two-components 
spinors. This way we keep control of the gauge couplings, which are explicit in the 
interaction vertices. This is justified since the couplings of the gauge eigenstates 
are dictated by the electroweak theory, but those of mass eigenstates depend on the 
amount of mixing [52]. Thus there are different, yet equivalent, bases of the mass 
eigenstates. We can choose a basis where the photino and the other gauginos reflect 
the couplings. Moreover, for our purposes it is enough to deflne just some generic 
mixing coefficients. In analogy with Eq. (l3.88|) . we can write 

X° = iV;.^°', z = l,---,4. (3.91) 

The non-diagonal mass matrix Y', that is analogous to (13.881) . can be indeed diago- 
nalized by a unitary matrix A^' so that N'*Y'N'~^ = Nj:,, with i,j running from 1 to 
4. This needs now to be rewritten in four components. The zino, the photino and 
the neutral higgsinos in the basis (I3.90p are respectively represented by the matrices 

-(t)' ^KS:) 

We define the four neutralinos as the Majorana fermions 

X° = (I ) • (3.93) 

By abuse of notation, the mixing factors for the spinors in Eq. (l3.92p will be labelled 
respectively by N~-, N'^^, Nj^. and A^^^, where the index i refers to the ith neutralino. 
The charginos are made up of gaugino and higgsino components, and mix in analogy 
with the neutral mass eigenstates. These can be defined in two-component formalism 
as [12] 

Xt = Vi.^l, xi = U^,^-, {1,3 = 1,2), (3.94) 



(3.92) 
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Figure 3.11: The diagrams which represent the amphtudes Mg and Mj. 



where i = 1,2. Expressing the W boson mass as 



(3.95) 



the charginos defined above are the eigenstates of the mass matrix (written in 2x2 
block form) 





^c-^ X 



with the 2x2 matrix 



X 



(3.96) 



(3.97) 



niw sin ^„ /i 
is written in the basis ip"^ = {W^, H^, W~ , H^), namely in components 

= (-zA+, = {-iX-, tPlJ, U = 1,2). (3.98) 



In this formalism, the mass terms in the Lagrangian have the form 

\- h.c. (3.99) 



and the mass eigenstates are related to the gauge eigenstates by the two unitary 
matrices U and V in Eq. (l3.94p : 



Xi 

xt 



V 



X2 



u 



(3.100) 



The mass matrix X is diagonalized as follows: 



WXV-^ = M, 



D 



^X2 



(3.101) 
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"a 



X 



+ 







' xj ' 


Ma 
r 





Figure 3.12: The diagrams corresponding to the amphtudes Mt and 



In four component notation, the charginos are the charged Dirac fermions 



Xi 



XI 
Xi 



X2 



Xi 

X2 



(3.102) 



w+ 



where by convention xi is heavier than X2- In the following, we will label with 
and ^Tjfy- the mixing elements of, respectively, the ith positive chargino with the 
wino and of the ith negative chargino with the wino W~ . 

After these preliminaries, we can now calculate the squared amplitude of the 
process f l3.85p . First consider the exchange of a gauge boson and the four-particle 
vertex, where Aj is the wave function of the i-th neutralino, as in Fig. (13. 11]) . By 
using Eg. (13. 91 p and the Feynman rules reported in the Appendix, we can write the 
amplitude of the diagram with exchange of TP as follows: 

, g cos Ow 



m 



AMp 
1 



nir 



[r/"^(A; - k'Y + r/^"(2A;' + kY - r/""(2A; + k'f] 



(3.103) 



It is evident from the above that in Mz the longitudinal modes in the TP propagator 
are cancelled by the SUGRA vertex, by the mechanism studied in Section I3.1.2[ 
Here we will not report explicitly the amplitudes of the other diagrams, which are 
shown in Figs. (13. 111) . (I3.12p and (I3.13p . Instead, we consider into detail the total 
squared amplitude M. of the process + W~ — )■ x° + G, which can be written 
as: 

M = M, + Mt + Mu + M, + Muiggs = Mstux + Mmggs . (3.104) 

We are interested in checking whether divergences such as (13.361) appear also in this 
case. The spin-summed square of the above depends on N^^^, Ofj, and it does not 
have a simple form. One can anyway fix some constraints on the mixing factors, 
since these are free parameters until they will be determined experimentally. First 
we report the squared amplitude of the four diagrams corresponding to Mstuxi then 
we consider also the Higgs bosons. By using the following equation. 



iVi-sin^w + A^Lcos^ 



W 



(3.105) 
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all the divergences which are proportional to l/m^m~ and to ^/fn'^ vanish. More- 
over, if 

2\N^f - 2N^^{0f; + Ofp + (|0,§p + lOgn = , U^^. = V^^, , (3.106) 

also those which are proportional to and to l/m^m~ are cancelled. If instead 

Ujyj^- and V^vj7+ are different, Eq. fl3.106p is rewritten as 

\Off + |Ogp = 2iV~ (O,^* + Oj* ) , U^^_ ^ V^^, , (3.107) 

and no anomalies are cancelled. In any case, the terms which factorize l/m~ still 
remain, 

0^ 1 1 
^-^ Mp cos^ 9w 3m ~ 

spin ^ G 



X 



\-\N^^ - Nt^smOwl'^j^ + Ai\N^f - N^^N^ Mw){s' + e + st)- 



-2N^^{0^* + Of*)s^ + 2NL sin 9ws^- 



-2 cos' Ow [2|iVp^J^t(. + t) + 2N^^[0\; + Of;)s' - AO^f^s'] | . 



(3.108) 



The above equation constitutes our net result. We considered Eq. fl3.106p for sim- 
plicity, and used = m-^y / cos 6y^. is the mixing factor of the neutral wino 
with the ith neutralino, in the basis given by 

{B,W\KM, (3.109) 

which is alternative to f l3.90p and contains the bino B. Both and B are linear 
combinations of 7 and 5, and the unitary matrix N diagonalizes the mass matrix 
r, see Eq.dMH). 

Consistently with our formalism, Eq.f l3.36|) is recovered in the limit when 9w — ?■ 
and Nyy-,Ojlj,0^ — )■ e'^*^, namely in the basis of SU(2)l gauge eigenstates. We 
remark that during this calculation no mass has been set to zero, and that no 
identities other than Eqs.f l3.105l) and fl3.106p have been adopted. 

The diagrams with exchange of Higgs bosons, which are reported in figure 
(I3.13p . contribute to the leading divergences only with the square, 

Esmn \Mm,,s\' ^ (/3 - «) + N'^. cos (/3 - a)p x 



X 



s' 2s' 



+ ^ r(2 cos^ Ow - 1) 



rri^im?- cos^ dwm'~ 



(3.110) 
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+ 





/\1 ' ' l^j 

Figure 3.13: The diagrams with exchange of the Higgs bosons 



where we have used the following formula for the Higgs masses [12]: 



m 



TTtr 



4m^m| cos' 



2/3 



(3.111) 



An evident analogy with Eqs. fl3.76p and (13.831) has been found. The scalars con- 
tribute to the factor l/m~ and reintroduce the divergences proportional to 1/m 



W 



and to l/m^m~. It can be shown that the general case U-yi^_ ^ ^jw+ "^o^s not 
change the result. This is evident, since (I3.110p depends only on s. Eqs. (l3.108|) and 
(13. nop complete the calculation in the basis of the mass eigenstates, and confirm 
what has been obtained in the basis of gauge eigenstates. 



It is instructive to check also whether the Ward identities (I3.17P hold in the 
present situation as well. In the basis of mass eigenstates, suitable constraints on 
the mixing factors eliminate the unwanted terms in the contractions of ka and k'^ 
with Mstux in (I3.104p . The scalar density S thus appears again. However, the Higgs 
bosons maintain their peculiar role, as it is clear from the following. 



k^M""^ = Sk'^ + Snk^ 
k'M^P = Sk^ + Suk"" 



(3.112) 



with the scalar amplitudes S and Sh written as 



S 



4Mf 



s — m 



(3.113) 



w 



and 



S 



H 



gmw 



[N^. sin {P-a){s-m]j) + N'g. cos (/3 - a){s - ml)] x 



y2Mp 

A,(l-75)V^^p''^ 



X 



[s-ml){s-mj^] 



(3.114) 



Violation of gauge invariance and therefore of unitarity is evident in Eq. (l3.11^ . by 
effect of Sh- Interestingly, the contraction of the longitudinal modes of the in 
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the propagator, that is 



s 



—I 



m 



2 

w . 



s - im\ 



V 



(3.115) 



with the 3-bosons vertex 



{rr\k 



k'y + + A;)" - r]'"'{2k + k'f] , 



(3.116) 



if non vanishing, would have provided with additional terms proportional to k'^ and 
k^ (or to k"' and /c'"). Thus there might have been a contribution to both the scalar 
densities S and Sh- These results are consistent with equation f l3.110p . 

To summarise, we have found that the squared amplitude of the scattering 



contains terms which violate the unitarity of the theory above a certain scale. The 
longitudinal polarizations of the Z° in the propagator vanish from the amplitude 
f l3.103p . and the Higgs bosons do not cancel the divergences. In this section, we 
have shown that it is not possible to restore unitarity by introducing the suitable 
scalar particles which are found in the observable sector of the MSSM. 

3.4 Conclusions and perspectives 

We have studied into detail the WW scattering in Supergravity in the broken phase, 
and compared our results with the case of massless gauge bosons. If in the high- 
energy limit the W bosons are considered massless a priori, as in Ref. |59]. the 
scattering amplitude does not contain any problematic terms. This is also the case 
of QCD p]. 

In contrast, the massive W provides with new structures, which could lead 
to the violation of unitarity in both the bases of gauge and mass eigenstates. In 
the annihilation diagram, the longitudinal degrees of freedom of the W boson in 
the propagator disappear from the amplitude by virtue of the supergravity vertex. 
This implies that the longitudinal polarizations of the on-shell W become strongly 
interacting at high energies. We show this effect in section I3.1.2[ where Eq. fl3.6ip is 
derived. We have also shown that in both the bases of gauge and mass eigenstates 
of the MSSM, no scalar particle can provide with cancellation of f l3.36p . 

Let us point out that in this chapter we do not claim that the theory of 
supergravity has inner inconsistencies. We have simply shown that when two massive 
bosons scatter and produce gravitinos, standard methods lead to contradiction of 
well-known theorems on the conservation of the supercurrent. 
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Cancellations of the molest terms would hopefully occur by adding further 
diagrams to the process. These can follow for instance from non trivial interaction 
terms in the supergravity Lagrangian (12.661) . which are not usually taken into ac- 
count in this kind of calculations. Moreover, the contribution of the hidden sector 
might account for the missing longitudinal polarizations. The latter possibility is 
the topic of current investigations. 

Our results are however phenomenologically interesting, if SUGRA is regarded 
as the effective limit of some unified theory. It is clear from Eg. (13.36]) that the SUSY 
breaking energy coincides precisely with the scale at which unitarity is violated. 
Accordingly, we can apply Eqs. (13.341) . (I3.108P and (I3.110p to particle physics and to 
cosmology. In fact, the process we have analysed can be observed at the LHC as a 
secondary reaction, for instance through gluon fusion. 

Moreover, we have shown in section (I3.1.2p that Eq. (l3.37p holds at any scale. 
Therefore our result contributes to the thermal gravitino production in the early 
universe, both at low and high reheating temperatures. How relevant the correction 
to the results of [53] and [59] would be, is the subject of a future work. 

We recall also that after inflation the MSSM flat directions, which are oscil- 
lating about the minimum of their potential, give a mass via their vev to the SM 
gauge fields which are produced by the infiaton decay p8|. Such gauge bosons scat- 
ter off each other in the primordial thermal bath, and recreate the phenomenology 
that is described in this chapter. Clearly, our results would affect this mechanism, 
providing with perhaps sizeable effects also in this particular cosmological context. 
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Chapter 4 



Scalar field oscillations and 
gravitinos 

In this chapter we discuss the cosmological evolution of a scalar field, which at the 
end of infiation begins coherent oscillations about the minimum of its potential and 
decays into gravitinos. By demanding that the gravitinos thus produced constitute 
the dominant component of Cold Dark Matter (CDM), bounds on the reheating 
temperature Tr are calculated. This is discussed in the paper [2]. 

We begin by considering the dynamics of the scalar field, that is here identified 
with the SUSY breaking field in the O'Raifeartaigh hidden sector of a model of gauge 
mediation with metastable vacuum. The scalar potential V{S), including radiative 
and supergravity corrections, is calculated. Studying the decay into gravitinos, we 
find that the reheating temperature can be easily as low as the BBN lower bound 
Tji ^ 1 MeV and still be consistent with gravitino Dark Matter from the decay of 
the GMSB hidden sector. 

4.1 Dynamics of a scalar in an expanding back- 
ground 

The cosmological evolution of the energy density associated with the oscillations 
of a scalar field was discussed long ago by Michael S. Turner and others [991 HOOj - 
Consider a scalar field (p with Lagrangian density 

£ = -^9^09>-\/(0), (4.1) 
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and a homogeneous and isotropic cosmology defined by the Friedman-Robertson- 
Walker metric Eq.f ll.ip . In such models the field is homogeneous as well, with 
equation of motion 

+ 3H0 + V'{(P) = , (4.2) 

where an overdot denotes d/dt and a prime d/dip. This is the equation of a damped 
harmonic oscillator, where the friction term 3H0 keeps track of the expansion of the 
Universe. Eq. (14. 2p can be recast as follows: 

^ = -3H(p + p), (4.4) 
d{pa^) = -pd{a^) . (4.5) 

As usual, a{t) is the cosmological scale factor. The energy density and the pressure 
of the field are respectively written as 

P = 1<P' + V{<P) , P = l^'- V{<P) . (4.6) 

The mass m of the field is initially much smaller than the Hubble rate. However, 
while the energy density of the infiaton drives the expansion of the Universe, H 
decreases with time. Accordingly, Eq.f l4.2p implies that when m ~ H, starts 
coherent oscillations about the minimum of the potential V{S), with frequency 
u ^ {4>/4>)- It is clear from the above that the energy density of the scalar field 
decreases, since the kinetic term (1/2)0^ is red-shifted away by the expansion. Let 
us assume that the frequency of the oscillations u is much larger than the Hubble 
rate. By defining 7 as the average of {p+p) over an oscillation and 7^ as the periodic 
part of 0^, we obtain 

0' = p + p= (7 + 7p)P- (4.7) 

Namely, the energy density of the scalar field decreases with time, since the oscilla- 
tions of are damped. This effect turns on particle production. 

The quantum particle creation can be studied by introducing suitable couplings 
which determine the total decay rate F of the scalar |100] . Following Turner |99] . 
we assume for simplicity that the particles produced by the decay are photons. Let 
be their energy density. We add a phenomenological term — 0^F to the rhs of 
Eq.f l4.3p . to obtain the evolution equations for p and p^, 

p=-(3H + F)7p, (4.8) 
p^ = -4Hp^ + 7Fp. (4.9) 
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Figure 4.1: Time-dependent shift Figure 4.2: The displacement from the 
of the minimum of the potential, minimum due to the time evolution of Sq. 



These can be easily solved to give the following: 



P = Po 

Pi ~ Pio 



ao 
a{t) 
ao 



e 

-4 



lT{t-to) 

a{t) 



+ Po 



ao 



-4 pVt 
JlVto 



ao 



-1 4-37 



''du . 



(4.10) 
(4.11) 



The subscript zero labels the value of that quantity at t = Iq, and u = ■yTt'. It is 
evident from the above equations that the energy density p during the oscillations 
decreases exponentially with time, and is transferred to the particles which are 
created through the deca}0. 

In this chapter, we study the production of gravitino Dark Matter by a mech- 
anism similar to the one sketched above. The scalar field, which we call S, is here 
included into a specific hidden sector in gauge mediation, and couples to the ob- 
servable sector via loops of the messengers [77]. A similar approach has been done 
recently in Ref. jlOlj . however our study of the dynamics of the S"- field is different. 
In |1U1] . the initial amplitude Sc of the scalar field oscillations is assumed to be 
determined by the Hubble corrections (normalized in units of the Planck mass) , 

AV ^ e^(3H2) ^ 3H2 

which imply \Sc\ = SfcM/lGvr^. This assumes that the inflaton dominates the energy 
density of the Universe when the S field starts oscillating. We will also include this 
contribution. However, in addition to fl4.12p . the scalar potential contains terms 
which are not H-dependent. As we will see, this implies that the minimum So{t) 



|5p-^(5 + St) + 



(4.12) 



If the particles created are nonrelativistic, all the 4's in Eas. (|4.10p and (|4.1ip become 3's. 
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is time-dependent in the presence of Hubble corrections. The S field tracks the 
minimum as it decreases towards the time-independent (H = 0) minimum. This 
generates a displacement of S from the minimum, as it is shown in Fig J4.1[ 

As H decreases, either the quadratic or the linear term becomes dominated by 
the H = factor. So the minimum of V is H-dependent until both \S\'^ and (5 + S^) 
become H-independent. The displacement from the minimum at this time, 6S, then 
determines the initial amplitude of the 5* oscillations about the time-independent 
minimum 5*0, as shown in Fig J4.2[ We first consider the dynamics of the GMSB 
sector generally, then we will apply it to the specific case of the model defined in 
Section 14. 2[ We consider a potential of the general form 

V = ^mlS^ -aS + - cS^ , (4.13) 

where the scalar field is assumed to take real values without loss of generality (with 
15*1 — S/\/2 for canonical normalization). When H decreases from a large value, 
either the linear or the quadratic term becomes dominated by the H-independent 
term. We define Hi and H2 by 

H = Hi when a = 3cH2 , (4.14) 
H = H2 when m| = Sli^ . (4.15) 

In the case H2 > Hi, the quadratic term becomes H-independent first. The potential 
when H2 > H > Hi is therefore 

V ^^mlS^ -ScR^S . (4.16) 

The minimum is then S'o(t) = 3cH^/m|. The field S will track this time-dependent 
minimum. If we define S = Sq + 6S, where 6S is a perturbation, the scalar field 
equation for S, 

dV 

S + mS = -—, (4.17) 

becomes the following. 



5S + m5S = -(Sq + 3H5o) - % 

oS 



JO 



5S. (4.18) 

So 



It can be shown that the higher-order terms in the expansion of the potential are 
negligible throughout our analysis. If the background is dominated by the inflaton, 
the Hubble parameter H depends on the scale factor a{t) as H oc a~^/^. Thus the 
above equation takes the form 



6S + 3H6S = -—^-ml5S, (4.19) 
2 rrig 
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and the displacement SS evolves as if it had a potential with minimum given by 

^=-V^- (4-20) 

This is illustrated in Fig J4.3[ 



V(S) 




Figure 4.3: Evolution of 58 in the potential with minimum at 882. 



On the other hand, if Hi > H2 then it is the linear term which becomes H- 
independent first. The potential when Hi > H > H2 is then 



V^-YC8^-a8. 
2 



(4.21) 



and the minimum is now at 6*0 (t) = a/3H^. In this case Eg. (14. 181) becomes 



which implies 



58 + 58 = ?,H^58, 

2 



3a 



58i 



2H2 



(4.22) 



(4.23) 



The shifts (I4.20p and (14.231) give the initial amplitude of the oscillations of 8 about 
the H-independent minimum in the two cases, where oscillations begin once H Hi 
and H H2 respectively. This modifies the gravitino density and the reheating 
temperature that were computed in |101j under the assumption that the initial 
amplitude is at H ^ ms- 

Before calculating constraints on according to the above dynamics, we 
have to include the field 5" in a specific particle physics model, in order to make 
phenomenological predictions. This is the topic of the next section. 
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4.2 Gauge- mediated supersymmetry breaking 

As we have seen in the previous chapters, the soft breaking terms in the matter 
sector must be of the order of the electroweak (EW) scale. On the other hand, 
the fundamental scale of spontaneous supersymmetry breaking (SSSB) can be much 
larger. Both with radiative corrections and at the tree level, the breaking occurs 
by new sets of fields and interactions, the "hidden sector". It is then transmitted 
to the observable sector in a variety of ways, depending on the interactions which 
communicate the breaking. For instance, these can be gravitational, corresponding 
to gravitational mediation. The SSSB scale is of order 10^^ GeV p ^lHlfl02] . In this 
case, the effective Lagrangian of supergravity contains non-renormalizable mixing 
terms between the hidden and observable sector, which are suppressed by powers of 
the (reduced) Planck mass Mp. Although this formalism is elegant and suggested 
by string theory, it is not clear whether it can solve certain problems such as flavour 
changing processes |1U2] . 

Another possibility is to mediate the breaking through gauge interactions, by 
introducing a set of new chiral supermultiplets, the messenger sector \103\ I104[ I105[ 
11061 1107| 1108] . Gauge Mediated Supersymmetry Breaking (GMSB) is particularly 
attractive, since it automatically avoids the problem of flavour changing processes. 
Moreover, it suits well to cosmology, since in this case the gravitino is generally 
the lightest supersymmetric particle (LSP), that can be Dark Matter in the case of 
conserved R-parity ^46j. The soft SUSY breaking terms are here introduced by the 
messengers. They couple directly to the hidden sector and indirectly to the matter 
sector through the SM interactions. In this case the gravitational interaction is 
automatically subdominant with respect to gauge couplings. 

If the model is renormalizable (like the MSSM), this theory is totally based 
on loops of the messengers. In the minimal model, these are a set of left-handed 
chiral supermultiplets q, q, i, i, which transform under SU(3)^ x SU(2)^ x U(l)y. 
They contain messenger quarks {ipg,ipg), squarks (g, g), messenger leptons {ip^.ipi) 
and sleptons (£, t) [H] . They must get very high masses in order not to have been 
already discovered, thus they should couple to a chiral singlet superfield S as follows: 

Wmess = hSU + k2Sqq. (4.24) 

Eq. (14.36 1 ) . that is discussed in the next section, has exactly this form. After that the 
scalar and the auxiliary components of S get a vev by supersymmetry breaking, the 
couplings above transmit this to the messengers, which become massive. It is not yet 
clear how the breaking occurs, however this can be achieved by either a dynamical 
SSB |109] or by virtue of an O'Raifeartaigh -type mechanism. The generation of the 
mass spectrum of the messengers is technically demanding and beyond the scope of 
this thesis. A detailed review can be found, for instance, in and in [41j. 
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Now, consider terms like (14 .24^ in the superpotential. After the SSSB, the 
scalar and fermion messengers get masses at tree level through mass terms which 
arise, respectively, in the Lagrangian and in the scalar potential. Loops of the 
messengers communicate this spectrum to the MSSM and generate the gaugino 
masses, as illustrated in Fig J4.4[ The q and q multiplets generate the gluino masses, 
while the loops of the i, i give masses to the wino and bino. The MSSM gaugino 
masses result as |109] 

(\ 

= -^A, (a =1,2, 3), (4.25) 

where the are the couplings of the respective interactions. The cutoff A, which 
fixes the normalization, is defined as 



A^{Fs)/{S). 



(4.26) 



The MSSM gauginos thus obtain their masses by virtue of spontaneous supersym- 
metry breaking, which has been communicated to the visible sector by one-loop 
processes of the messengers. 

In contrast, scalars do not get radiative corrections to their masses from one- 
loop diagrams, but from processes at two-loops. These are much more complicated 
than those in Fig J4.4l and involve messenger fermions and scalars, together with 
ordinary gauge bosons and gauginos. The mass of each MSSM scalar results as 



rrij 



2K' 



^^^^ 



1,2,3) 



(4.27) 



where the Ca{i) are the quadratic Casimir invariants [16]. Eqs. (l4.25|) and (I4.27P 
are valid in the limit of mass splittings within each messenger multiplet which are 
small with respect to the overall messenger mass scale, that is the case of the most 
common models. 
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4.2.1 Gauge mediation in metastable vacua 

We now introduce a particular class of GMSB models, which have been recently 
introduced by Murayama and Nomura in [1101 1111] . 

Such models are defined by three basic ingredients: i) the Kahler potential 
contains a negative quartic term in the supersymmetry breaking field, ii) an acci- 
dental U(1)r symmetry in the superpotential terms which are determinant for SUSY 
breaking, iii) an explicit messenger mass term. One can see this by writing down 
the superpotential |111] . 

W = -^'S + kSff + Mff, (4.28) 

with no hidden sector yet specifiecil. 5* is a gauge singlet chiral field, that can be 
either elementary or composite. The parameter k can be taken real and positive 
without loss of generality, the scales n and M have dimension of mass, and the mes- 
sengers / and / belong to 5 + 5* representations of SU(5) to include Standard Model 
interactions. The Kahler potential for the messengers is assumed to be canonical, 
whilst for the chiral field S it looks as follows, 

1014 

K=\S\'-^-i. + On^] , (4.29) 



4A2 \A\ 

after expansion around 5 = 0. A is a mass scale which has to be assumed as a 
physical cutoff for the low-energy theory, in the absence of an explicit hidden sector. 
An approximate low-energy U(1)r symmetry on S makes it possible to obtain this 
form of the Kahler potential (this will be shown later in our particular model). For 
instance, the first two terms in (14.281) may have such symmetry if S, f and / carry 
respectively the charges 2, and 0. The third term in (I4.28P violates this symmetry 
anyway, but we will see that its effect on the Kahler potential can be neglected under 
certain assumptions. At the tree level the potential is given by 



v=W-kff\' 



A2 V A4 
and has a global SUSY minimum at 



+ \kS + M\\\f\' + m, (4.30) 



S = -^, f = f = ^. (4.31) 

There is also a local SUSY breaking minimum (which is also a metastable vacuum) 
atS' = / = / = 0, if > kjji^. The latter condition is also necessary to avoid 



^An O'Raifeartaigh hidden sector will be included in the next section. 
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tachyonic messengers, as the scalar components of the fields have masses m| = /i^/A^ 
and m'j = ± fc/i^. Since, as found in the previous section, the messenger loops 
generate the gaugino masses as 

" 167r2 M ' ^ ^ 

by demanding these to be of order 100 GeV to 1 TeV we find that 

4^ ^100 TeV. (4.33) 
Moreover, the gravitino mass is given by (see Chapter [2]), 

{Fs) . X 

^ Mp Mp' ^ ^ 

therefore we find that if gravity mediation is subdominant (mg ^ 10 GeV), the 
SUSY breaking scale /i is smaller than 10^° GeV. 

A detailed discussion about the U(1)r symmetry in the superpotential f l4.28p 
can be found in the paper by Murayama and Nomura Ref. [lllj . Let us anyway 
sketch here some remarks. The SUSY breaking minimum at S = f = f = might 
be viewed as a consequence of this symmetry, so it is clear how it plays a central 
role in the class of GMSB models we are considering. Nevertheless, violation of 
such symmetry may come from several effects. First of all, from nonlinear terms 
in S, such as S"^ or S^, that may appear in the superpotential. In the case of an 
elementary S field, one must simply assume that they are somehow negligible. On 
the other hand, if S is composite, these terms can be automatically suppressed. 

Another possibility of violation comes from loops of the messengers, which 
violate U(1)r by virtue of the mass term in the superpotential. This induces a 
linear term in the Kahler potential, which makes the SUSY breaking field deviate 
from the origin. As we will see in the following, such effect is responsible for the 
decay of the field S into gravitinos. 

We remark that, at this point, is it quite evident that the models of gauge 
mediation with metastable vacua so constructed have only four parameters: fi, k, A, 
M. Remarkably, the gravitino mass range allowed by this framework is very large, 
1 eV ^ 10 GeV. We will go further in the analysis of this formalism in the 

next section, where a particular model with an explicit hidden sector is addressed. 



4.2.2 Model with an O'Raifeartaigh hidden sector 

As we stated in the previous discussion, the form of the Kahler potential Eq.f l4.29]) is 
generated by the radiative corrections of the fields in the hidden sector. Moreover, in 
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order to have a phenomenologically viable model, the cutoff A should be expressed 
in terms of the physical parameters of such hidden sector, which then needs to be 
introduced explicitly. Let us then consider a field content of the O'Raifeartaigh 
type, with superpotential 

Wh = -n^S + A^X^ + mXY , (4.35) 

where the singlet fields S*, X and Y have a canonical Kahler potential at tree level 
(up to terms suppressed by the Planck mass). Here S is elementary, thus terms like 
S"^ , , , SXY , SY'^ are assumed to be negligible . By recalling the 
general superpotential introduced in Eq. fl4.28p . we then have: 

W = -fi'^S + XSX^ + mXY + kSff + Mff = Wh + Wm. (4.36) 

Along the direction of the field space corresponding to X = 0, F = 0, we still find 
a global SUSY minimum for the scalar potential at 

and a local SUSY breaking minimum (which is a metastable vacuum) at 5* = / = 
/ = 0, if > kfi"^. In fact, since at this point there is incompatibility between 
Fs = and Fy = 0, the field S breaks supersymmetry. We have seen in section [2. 1.21 
that in O'Raifeartaigh models the SUSY breaking field S is usually a fiat direction 
of the potential at the tree level. In fact, the classical scalar potential of the hidden 
sector, that is obtained from the F-terms of the superfield multiplet, holds as follows 

V,\S) = |Fs|' + |FxP + |Fy|2 = / + A2|X|^-A/i2(x2 + Xt2) + 

W(|X|2 + |r|2) + 2Xm{SXY^ + S^X^Y) + AX^\X\^\S\^ . (4.38) 

Thus in the direction X = Y = 0, the above is constant (or zero, by rescaling) for 
any value of 5*. However, the radiative corrections originated by loops of the X, Y 
fields and of the messenger fermions and scalars, have the property to lift this flat 
direction. We discuss this effect in the next section. 



4.2.3 Coleman- Weinberg potentials 

In 1973, Sidney Coleman and Erik Weinberg published a paper, where they claimed 
that radiative corrections may produce spontaneous symmetry breaking in theories 
in which such breakdown is absent at the tree level [112] . Clearly, they could not 
consider supersymmetric theories at that time. However, starting from a toy model 
concerning massless scalar electrodynamics, they have shown that the effect is valid 
for both abelian and non-abelian gauge field theories, if the couplings are small. 
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The extension of this effect to the Standard Model and to supersymmetric the- 
ories is naturally achieved. In several fields of research (i.e. in solid state physics, or 
in the study of phase transitions), this result has proven to be extremely important. 
In particle cosmology, there have been applications in many fields, and in particular 
in theories of inflation [131 111 1113] . 

In this section, we show that the Coleman- Weinberg effective potentials for 
our model, which are generated by loops of the hidden sector and of the messengers, 
produce a series of important effects for the phenomenology. Namely, they i) lift the 
flat direction in the classical potential V{S) by generating a mass term for the field 
S; ii) provide with an explicit form of the otherwise arbitrary cutoff A in terms of 
the parameters of the GMSB model; iii) induce a linear term in the Kahler potential, 
which through supergravity corrections enters V{S) and makes S deviate from the 
origin. 



Contributions of global supersymmetry 



The mass matrix in the basis of the X and Y fermions ipx and ipy is given by 

(4.39) 



Mi 



2XS m 
m 



while the F-terms of the hidden sector, in the basis {X, X^, Y, F^), provide with the 
following mass matrix for the scalars |111] : 

/m2 + 4A2|^|2 -2A/i2 2XmS^ \ 

-2Xfi^ m^ + AX^\S\^ 2XmS 

2XmS 
2XmS^ 



V 



m 



(4.40) 



Given the mass matrices and by equations ( I4.39P and f l4.40p respectively, the 
radiative corrections to the scalar potential can be now calculated with the following 
formula 

1 



eff 



647r2 



TrfMjlnMcj] 



for the scalars, and 



Tr[(|Af^|^)^ln(|M^n] 



(4.41) 



(4.42) 



for the fermions |112] . The minus sign comes from the Fermi-Dirac statistics. The 
formulas (I4.4ip and (I4.42p give the following Coleman- Weinberg potential for the 
hidden sector. 



Veff{S) = Vff^{S) + 2Vf^f{S) 



4:, A 



X^li 
37r2m2 



|5p 



6, ,4 



3AV 
107r%n4 



\st + ... 



(4.43) 
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where we have dropped a constant, and kept the dominant terms in A/i^/m^ and in 
S. The above equation agrees with Since this contribution has to be added 

to the tree level potential in Eg. (14.381) . it has been verified that at the quantum 
level, the field S is stabilized at the origin of field space (because of the expansion 
in S). Moreover, we see in (I4.43P a positive mass term, which for small values of S 
dominates over the quartic term. This is interesting from the point of view of the 
dynamics of the model. 

Also the one-loop corrections to the Kahler potential are important for phe- 
nomenology. They can be calculated with the following formula 



7^(1) 



327r2 



:Tr 



In- 



VI 



A2 



(4.44) 



that is a particular case of the more general expression 
1 



327r2 



rTr 



327r2 



/^/iln^-2Mln- 
Tr Mi ( In -^-1 



2Tr 



In 



Ml 
A2 



- 1 



(4.45) 



which is valid for any renormalizable supersymmetric theory with = 1 |114j . Here 
M| and are the mass matrices of, respectively, the chiral and vector superfield 
sectors. In particular. 



for the chiral sector and 

Mab 



or 



(4.46) 



(4.47) 



for the vector sector (which is however absent in the model here discussed). In the 
case of messenger loops, we multiply the right hand side of Eq.f l4.44]) by a factor of 
5, to take into account for the representations of SU(5). 

In the previous section we have found that, in this class of GMSB models with 
metastable vacua, as long as one does not introduce an explicit hidden sector, the 
non-minimal Kahler potential for S is written as: 



K=\S\' 



\S\^ 
4A2 



+ 



A4 



(4.48) 



In the case at hand we have defined the hidden sector, therefore the physical cutoff 
A can be written in function of the parameters of the O'Raifeartaigh model. By 
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using Eg. (14 .44 p . for small S and at the lowest order in A/i^/m^, they result as 



87r2 



151 



A^ 



127r^m^ 



151 



which correspond to the potential 

K{S) ^ \S\^ 

with a cutoff 



\si 

4A2 



+ 



(4.49) 



(4.50) 



' (4-51) 

since A^/Svr^ ^ 1 even for a strong coupling A ~ C^(l)- Remarkably, quantum 
corrections have generated both a mass and a non-canonical Kahler potential for 
the supersymmetry breaking field. 

The contribution of the messengers is fundamental to our discussion, as they 
generate in the potential linear terms in S. These are responsible for the displace- 
ment from the origin that induces the oscillations of the field S, which is essential for 
cosmology and will be considered in the next sections. The scalar and fermion mass 
matrices and M^, in the bases [f^P^J^P) and {jpf^tpj), can be respectively 
written as 

/ |M + A;5p 

\M + kS\'^ -kij,^ + kXX^ 

-kfi^ + kXX^^ \M + kS\^ 

\ -kfi'^ + kXX'^ 







-kfi^ + kXX^^ \ 



\M + kS\^ 



and as 



M + kS 
M + kS 



(4.52) 



(4.53) 



At the lowest order in k^'^/M'^ and for small values of S, the above matrices give 



5M2 
16^ 



(4.54) 



that agrees with the results of [lOlJ (notice the linear term in S). The F-terms of 
the messengers give the following tree-level potential: 



(4.55) 
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The Coleman- Weinberg potential of the messengers is obtained by summing the 
contributions of the loops of the messenger scalars and fermions, as before. Keeping 
only the dominant terms in fc/x^/M^, in A/i^/m^ and in S, we find 



X 



16772 

^247r2M6 



{S + S^) 



_5kY_ 
327r2M2 



{S^ + 5t2) 



+ ^(X^ + + 4|X|^) + Xfi\X' + Xt2) 



151 



(4.56) 



where we have dropped an unimportant constant. This expression agrees with the 
literature |10H 1111] for X = Y = 0. Note that the messenger loops generate a mass 
term for the field S as well. 



The radiative corrections of global SUSY to the scalar potential, i.e. Eqs. fl4.43p 
and f l4.56p . which are generated by loops of the supermultiplets in both the hidden 
and messenger sectors, can be now added to the tree-level potentials (14.381) and 
(I4.55p . Along the direction X = y = / = / = Oin the field space, where both the 
SUSY and the metastable minimum are localized and where the 5* field oscillates, 
the classical plus quantum potential can be rewritten as follows. 



Vsusy{S) 



4 , m 



2 3A/i^ 15A;2/i^ 5A;2/i4 
+ + o + —V In 



327r2 



+ 



167r 
5A;V 
'327r2M7 



167r2 R^ 

5A;V^ 
487r2M6 



151 



(5 + 5t)|5p 



3AV^ 15A;V 



+ 



107r2m4 647r2M8 



\s[- 



(4.57) 



where i? is a renormalization scale. This is added to the supergravity corrections 
which are discussed in the next section. We then obtain the full scalar potential 
V{S) that is used to calculate the gravitino abundance. 

Supergravity corrections to the potential 



To embed the model of gauge mediation in supergravity, we can use the general 
formula [35] 



VsuGRA = e''p 



K 



E 



d(t>ad(t)l3 



\wl_ 

Ml 



where Mp = 2.43 x 10 GeV is the reduced Planck mass and 



DW 



dW W dK 



p ^'Ta 



(4.58) 



(4.59) 
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is the Kahler derivative of the superpotential W. At this point, it might be inter- 
esting to show first the effect of the quartic term in 5* in the Kahler potential at 
the tree level. We will then include the radiative corrections from the loops of the 
messengers and of the hidden sector. Eg. (14.581) for the SUSY breaking field is here 
rewritten as 



V, 



SUGRA 



1 



\si 

A2 



2\ -1 



X 



Ml 



Ml 



Ml 



(4.60) 



with the Kahler potential fl4.50p . At the lowest order in l^p, this is expanded as 

12 |^|4 



ysUGRA{S) ^ /X 



A2 



4A2M|, 



/ ( 1 + — ^ + 



127r2 m^Ml 



(4.61) 

The term proportional to 15*^ has a positive coefficient, consistently with the liter- 
ature |101j . Using the radiative corrections (I4.49P and (I4.54p . expansion in S gives 



K^^\S) = K^^\S)+kI^\S) ^ \S 



^Af(5+St)-J^^(S+5t)-^^ 



167r2" " "~ ■ ~ ' 327r2 M 
at the lowest order. Thus the following dominant terms are generated. 



127r2 

(4.62) 



Vsugra{S) ^ /i- 



\S\' f25P 



P I 647r4 



M 

wi 

'mi 




+ . 



(4.63) 



where we write only the terms which are linear and quadratic in S. This agrees with 
f l4.6ip in the limit /c — )■ 0. We can further assume that 



25A;2 
64^ 



M 
Wp 



37r2 



Mr 



m 



(4.64) 



therefore the supergravity corrections can be rewritten as 

Vsvgra{S) ^ /.^ \^,\S? + -|^(5 + 5t) - 2^ 



+ ... 



(4.65) 



Clearly, the sign of the term that is linear in S is now determined by the coupling 
constant k. In particular, a form of the potential such as 



Vsugra{S) ^ fi' [+c'\S\' - d\S + S^) + ...], 



(4.66) 
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would be recovered if the following condition is satisfied: 



k <V2 



M 



m 



5.77 

> 

ms as 



G 



(4.67) 



The above equation directly relates the gravitino mass to the gluino mass 
and to the strong coupling constant a^. It is calculated by using Eq.f l4.74p in Section 
14.31 However, fl4.67p is not consistent with gauge mediation, where the gravitino is 
the LSP. In particular, the framework that is used here admits gravitino masses 
10 GeV, thus we conclude that also the coefficient of the linear term in ( I4.65P 
is identically positive. 

The supergravity corrections to the scalar potential of the model can now be 
written by adding ( I4.65P to the potential for S which is generated by the infiaton, 
whose energy density dominates the Universe before the field starts oscillating |101j . 



ViS) 



3H' 



hk 



M{S + S^) + 



(4.68) 



This provides, at the lowest order in the SUSY breaking field, with 

5k 



Vsugra{S) 



3W 



151 



167?^ 



M{S + S^) 



+ 



. (4.69) 



The above supergravity corrections, together with Eq.f l4.57p . are used in the next 
section to calculate the gravitino abundance and the reheating temperature. 



4.3 Gravitino Dark Matter and reheating tem- 
perature 

The results of Sections |4. II and 1^12) can be now applied to gravitino production. Here 
we calculate the gravitino number density as a function of the reheating temperature 
Tji and of the parameters of the GMSB model. We find bounds on Tr by demanding 
that the gravitinos produced by the decay of the scalar field constitute the dominant 
component of Cold Dark Matter. The scalar potential of the field S is obtained 
from the superpotential fl4.36p and from the Kahler potential f l4.29p . as a sum of a 
classical part and the corrections generated by quantum loops and by supergravity. 
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V{S) = VsuSY + VsUGRA^{SR' + ^,+-^l^]\S\' + 



By summing equations (I4.57P and fl4.69p . we obtain 

3^m2 ' 487r2M6 
+l|.M(S + «')(-3H'-f + ^). (4.70) 

where we have kept only the hnear and quadratic terms in S, in the direction 
corresponding toX = y = / = / = 0. These are dominant for sufficiently small 
values of the field. 

The scalar field can be treated as a modulus, with interactions given in [77] . 
The decay into gravitinos is not helicity suppressed, and gravitino production from 
5* decay dominates production from inflaton decay. The mass scales here considered 
are indeed consistent with the discussion in Ref . |68] . 

Dominant decay into a pair of gravitinos {S — ip^'^u) occurs if the cutoff A 
satisfies [101] 

A<8xlO-,/^f^)"V^) GeV. (4.71) 
Vo.lVlTeV/ VI GeV/ ^ ' 

The inclusion of the O'Raifeartaigh sector makes it possible to relate this directly 
to the mass scale m in the hidden sector using Eq.f l4.5ip . 

™,1.5xlO»A^GeVxy|I(^)-'(^). (4.2) 

This corresponds to an upper bound in the range 10*^ GeV — 10^^ GeV for a gravitino 
mass between IkeV and 10 GeV (the lower bound on is discussed in Sect J4.4p . 
The messenger mass scale M is related to the gluino mass and to the SUSY 
breaking scale fi by [lOlj 

M = ^ , (4.73) 
therefore, by using Eq. fl4.112p . we find 

This implies 10^ GeV ^ M ^ 10^*^ GeV for the gravitino mass range reported above. 
Together with (14. 72 p . it also gives the following approximation for the coefficient of 
the linear term in the scalar potential fl4.70p . 



2^^ 2P^n ^ / , 2fcV' , 
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over the entire gravitino mass range considered. Similarly, for the coefficient of the 
quadratic term, 

We now briefly comment on the thermalization of the hidden sector, which was 
discussed for example in [116j . For Tr smaller than the messenger mass M, the 
scattering of S particles by thermal MSSM particles will be due to messenger 
loops, with rate ^ a^T^/M*^ (a^ accounts for couplings constants and numer- 
ical factors). This can be compared with the expansion rate at T < Tr, namely 
H ^ T^/Mpi. Accordingly, the hidden sector will not be thermalized as long as 
Tr < {M/asMpiY^^M. Then f l4.74p implies that the corresponding upper bound 
on the reheating temperature for which the S field is not thermalized is 1 TeV to 
10^2 GeV (using aj^ ~ 1). 

One should also consider what happens when adding a constant term ~ m^Mp 
to the superpotential of the model. This is needed to tune the vacuum energy to 
zero in supergravity [351 HZ]- Recalling Eq. fl4.36p . we can write [77] 

W = -iJ,^S + mQMj, + ... (4.77) 

This generates both linear and quadratic terms in S through the factor |PFp in the 
supergravity corrections fl4.58p . that is 

VsuGRA = eMK/Ml) (^.. - 3^^ ^ + Q^J?m^{S + ^t) 



M|, 167r2 M|, 

,4 



My ' 4v^7r2 Ml ' ' 



(4.78) 



Eqs. (14.751) and (14.761) change accordingly: 



2A;V 32V37rV^ 1 

5^MMp ' ' ^^-^^^ 

3H^ + TT-A--^T7T • (4-80) 



37r2m2 4v/37r2 M|, 

For the entire gravitino mass range here considered, the contribution of the addi- 
tional term in f l4.80p is always negligible. On the other hand, if k ^ 0.1 the new 
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term in (14.791) dominates for certain values of mg, thus we take it into account as 
well. The scalar potential of interest is therefore the following, 



37r2m2 J 2 Sn'^ \ M"^ 5k MMp 

(4.81) 

where it is assumed that S is real for simplicity, which will be the case on minimizing 
the potential if all parameters are considered real. Equation (14.8 ip can thus be 
rewritten as follows. 

Comparison of this expression with the potential 

V = ^mlS^ -aS + 3H2 - cS^ , (4.83) 

gives the parameters which determine the dynamics of the field S, as discussed in 
Section 14.11 It is now possible to calculate the gravitino abundance from the decay 
of the scalar field S. 

The oscillations start when the Universe is still dominated by the inflaton, thus 
T > Tr. The energy density can be written as 

P = g^^T' , (4.84) 

where T is the photon temperature and g^, is the number of effectively massless 
degrees of freedom. If the spacial curvature is /c = 0, the Hubble rate in radiation 
domination can be written as 



H(T) = W^^. (4.85) 
^ ^ V 90 Mp ^ ^ 

The entropy density of the Universe is dominated by relativistic particles, therefore 
p = p/3 and (14.841) gives 

47r2 

s = g^s T^- (4.86) 

^ 90 ^ ' 

Assuming thermal equilibrium, the conservation of entropy per comoving volume, 
namely s{T)a^(T) = constant, can be rewritten as 

g^{T)T'^a^{T) = constant , (4.87) 
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where a{T) is the scale factor at temperature T. Let ns osc be the number density 
of the scalar and aosc the scale factor at that time. The number density ns{T) at a 
generic temperature is given by 



3 / \ 3 / /m \\ 3 



B.{Tr)V g(T)T^ 
^ ^ ^ ^ nsosc. (4.88) 



since H oc a~^/^ during inflaton domination. Recalling Eq.f l4.85|) . the above becomes 

(rr\ T^"^ 9{T)T^TRnsosc , . 

Therefore the number density to entropy density of 5* is just 

ni^jR_nsosc . . 

s AMI, HL • ^ ^ ^ 

From the scalar potential Eq. fl4.82p . ps = m|S'^/2, which implies ns osc = Ps/^s = 
{l/2)msSl^^, where Sosc is the value of the field at the beginning of the oscillations. 
Therefore 

s 8M|HL • ^ ^ 

Since each 5* scalar decays into a pair of gravitinos, the gravitino number density to 
entropy ratio is given by 

ng _ msTRSl 



'osc 



(4.92) 



where the initial oscillation amplitude around the minimum, Sosc, is given by SSi or 
SS2 at the time when all H dependence becomes negligible. 

Consider first the case where H2 > Hi, namely where H^^,^ ^ B.^ = —a /3c and 
Sosc = SS2. Therefore 

n^_msTRlS^ 

s %Ml H? • ^ ^ ^ 

The parameters a and c and the mass nis are fixed by the messenger and hidden 
sectors. 



4v^7r2M y Mp 
5kM 
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Using these we obtain 



20257r^ k^mjm^TR 



I6V2 X^^airrvLMl [ 



127ra3 / mg 



(4.97) 



where we have used fl4.73p and fl4.112p to ehminate M and /x. Therefore 



^ - 0.11 X ^ 
s Vo.l 



^3 y [ \ 

ITeVJ \lGeVJ 



m 



1014 GeV 



108 GeV 



1 + 



127ra 



3 / 



5P 



(4.98) 



The (rescaled) density parameter of the gravitino Og/i^ (see Chapter [T]) can be now 
calculated as follows. The entropy density is s = 27T'^g{T)T^/45, with at present 
g{T) = 2 (for the photon). Since the CMB temperature is T = 2.4 x lO'^^ GeV, 
and the energy density of the gravitino is = mQ{n^/s)s, we obtain 



(4.99) 



with the present critical density pc = 8.1 x 10 ^'^/i^ GeV^. Using the above formula 
we find 



VI GeV/ 



m 



1014 GeV/ V 108 GeV 



R 



1 + 



127ra 



5F 



3 / ^G 



(4.100) 

By demanding that the gravitino is the principal constituent of Cold Dark Matter, 

fig/i^ ^ C(O.l), (4.101) 

we then obtain the reheating temperature Tr as a function of the parameters of the 
model 



A 



14 



Tr ^ 0.6 GeV X — — 
^ Vo.l 



"3 



^3 y 
TeVy 



m 



1014 GeV 



5A;^ \ mg 



IGeV/ 
1 -3 



(4.102) 



In the case Hi > H2, H^^^ ^ Rl = m|/3 and Sosc = ^Si. We then obtain 



ns 



182257r3 /t^m^m^Tjj 
256v^ AiOaim|M|, f ' 5^ 



(4.103) 
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Therefore 



and accordingly 



^ 3.9 X 10"^— ( — 



X 



m 



1014 GeV 



"3 



Aio VO.l 



-2 



1 Tev) 



108 GeV / 



1 + 



IGeV 
127ra3 



n 2 



(4.104) 



5.9 X 10^ 



A 



10 



X 



m 



Q3 

0.1 



lTeV7 



1014 GeV/ V 108 GeV 



IGeV/ 

127ra;3 



1 + ^7^^ 



5k 



ma 



(4.105) 



The reheating temperature Tr required for the correct density of gravitino Dark 
Matter is then 



T 



R 



ITGeVx^p 
k^ Vo.l 



m 



1014 GeV 



1 + 



ITeVJ 
127ra3 
5A;2 



-2 



■V 



IGeV/ 

-2 



(4.106) 



Eqs. fl4.102p and (I4.106P are our main results. We see that the model can accom- 
modate a very wide range of Tr and still be consistent with gravitino Dark Matter 
from the GMSB sector. The most striking feature of these results is their extreme 
sensitivity to the parameters of the model, in particular the O'Raifeartaigh sector 
coupling A. This means that the gravitino density from decay of the SUSY breaking 
scalar in the GMSB sector can account for Dark Matter for essentially any value of 
the reheating temperature above the BEN bound, Tr >, 1 MeV. 

Comparing the above with the previous results of [lOi j , which assumed that 
oscillations about the minimum began at Hose ~ t^s with Sosc ~ V^c (for real 5"), 
we find 
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where ng/s|„„ is the value given in |101] . Namely, 



(4.107) 
(4.108) 



^ 2 X 10-'%^ ( ^ 

ET VO.l 



■V 
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^3 

ITeVy VlGeVy V lO^ GeV 



A 



1014 GeV 



(4.109) 
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If Hi > H2, |a/cm|| > 1 and so there is a strong enhancement. Similarly, if H2 > Hi 
then there is a strong suppression of the gravitino abundance. Note that while in 
the existing literature Tp> is necessarily written in terms of an arbitrary cutoff [lOlj . 
in our analysis there is a fully defined perturbative hidden sector, eliminating the 
cut-off in favour of the masses and couplings of the model. Therefore what could 
have been a limitation, namely focusing on a specific O'Raifeartaigh -type sector, 
is in fact an advantage. Since there is no cut-off, there is no need to constrain 
the parameters of the model in order to make it consistent with the validity of the 
effective theory. 



4.4 Constraints on the GMSB model 

Here we derive various bounds on our model, first by demanding consistency of the 
series expansion of the Coleman- Weinberg potential. The coupling constants A and 
k and the mass scale m are indeed constrained as follows. The scalar potential 
Eg. (14. 70 p is obtained by expanding in terms of A/x^/m^, so the condition 

^«1, (4.110) 

must be satisfied. Therefore 

m>A^/V- (4-111) 

Now we derive further constraints on the mass scale m of the hidden sector, and show 
that the gravitino mass range which is consistent with the model can be accordingly 
constrained. First, by using the relation between and mg, 

= y^m^Mp , (4.112) 
the lower bound f l4.11ip can be recast as 



m > A^/V = 2 X lO'X^' GeV^^^ . (4.113) 

In Section 14. 3[ we assume that the gravitino is the main decay product of the field 
S when this rapidly oscillates around the minimum. Accordingly, the upper limit 
Eq. flTO]) holds: 



.5 X 10^^ A^ GeV x (^) (^) . (4.114) 

V 0.1 VlTeV/ VlGeV/ ^ ' 



Clearly, fl4.113p and f l4.114p must be simultaneously satisfied. This means that the 
formulas f l4.102p and fl4.106|) for the reheating temperature are valid only for certain 
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gravitino mass ranges, depending on the values of A. In particular, 

A ^ 10-2 ^ IMeV < mg < lOGeV, (4.115) 
A ^ 10-^ ^ 1 keV < mg < 10 GeV , (4.116) 

where we have assumed only that ^ 1 TeV. Remarkably, since these results follow 
from the consistency of the Coleman- Weinberg potentials and from phenomenology, 
they do not rely on any approximations. 

Next, let us calculate constraints on the coupling k and show that it should 
not be too small. In analogy with Eq.f l4.1lTl) . the effective potential fl4.56p for the 
messenger sector is valid if 

(4.117) 

By recalling Eq. fl4.74p . i.e. 

M«3,3xlO'=GeVx.(|i)(i^)(^). ,4.118) 
we can substitute the above equation and ( 14.1121) in (14.1 17p . Thus we finally obtain 

lO^^/aiV^^y.^x-. 

2k \aj WtcV/ VlGeV/' ' ' 

This condition is satisfied for k ^ 0.01 in the entire gravitino mass range in equations 
( 14TT511 and (HTTHll . 

One should also take into account the free streaming length Xps, namely the 
distance that the gravitinos produced through the decay can travel until they become 
non-relativistic. If we want them to be Cold Dark Matter (CDM), this quantity must 
be smaller than (9(100 kpc) [115] . With the notations of Ref. |101j , the comoving free 
streaming length at matter-radiation equality can be defined as 

v~(t) 

\ps= / -^dt, (4.120) 
Jtn «(^) 

where a{t) is the scale factor, is the time at the S decay and teq ~ 2 x 10^^ sec is 
the time at matter-radiation decoupling. The velocity of the gravitino is given by 

ms f ao 
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where is the scale factor when S decays, and we have assumed nis ^ mg. 
Integrating fl4.120p and recalhng that = Svr^m^/ A^, one finds 



3 

2 I , _ , 3 . ,3 



where g^: is the number of relativistic degrees of freedom at the time of decay. Since 
the free streaming length must be smaller than 100 kpc, one obtains a lower limit 
on the gravitino mass, 

m. ^ 5.44 X HZ f A^)^ f,^) GeV. (4.123) 



■ A2 VlOO/ VlOi^GeV 

Clearly, depending on the coupling A, the above equation determines a set of lower 
bounds on the gravitino mass. For A = 10~^, the former constraint ^ 5MeV 
still holds, whereas for A = 10^^ we find a slightly different mass range. In this case, 
the gravitinos produced by the S decay are CDM only if 

A ^ 10-^ 22 keV < mg < 10 GeV. (4.124) 

By taking a look at Eg. (14. 1231) . we see that this is not the most stringent mass 
bound. The lower limit ^ 22 keV is indeed a function of m, which in principle 
can assume higher values, as soon as Eq. (l4.114p is satisfied. 



4.5 Conclusions and outlook 

In this chapter we have considered the dynamics of the SUSY breaking scalar S* in a 
GMSB scenario with metastable vacua and the production of gravitino Dark Matter 
through its decay. Our results for the cosmological evolution of the scalar field and 
the resulting gravitino density are significantly different from previous investigations. 

We have shown that since the Universe is expanding, the minimum of the po- 
tential 5*0 is time- dependent. The 5* field tracks the minimum, which generates a 
displacement 6S. Once the minimum becomes H independent, 5* begins coherent 
oscillations about the minimum with initial amplitude determined by the displace- 
ment. This produces a very different S oscillation amplitude and so gravitino density 
as compared with previous analyses jlOlj . By considering a generic potential V{S), 
we have shown that there are two possible values of SS, depending on whether the 
quadratic or the linear term first becomes H-independent. The resulting gravitino 
density can be highly suppressed or enhanced as compared with the previous esti- 
mate of [101] . A striking feature of the gravitino density is its extreme sensitivity 
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to the parameters of the model, with the reheating temperature having a or A^^ 
dependence on the superpotential couphng of the O'Raifeartaigh sector. It is there- 
fore easy to account for gravitino Dark Matter with an arbitrarily low reheating 
temperature. This could be significant for the cosmology of GMSB models, since a 
low Tji is consistent with several cosmological scenarios |117j . 

One might also be interested in completing the model with a source of baryo- 
genesis. A natural possibility would be Affleck-Dine baryogenesis. Indeed, for 
~ 1 GeV, Affleck-Dine baryogenesis combined with Q-ball decay in GMSB could 
account for both gravitino Dark Matter and the baryon asymmetry |118] . but for 
smaller mg it would provide only the baryon asymmetry. As an additional remark, 
we note that the O'Raifeartaigh GMSB sector is a well-defined perturbative model. 
Therefore the decay of the SUSY-breaking scalar field is addressed without cutoff, 
nor with any stringent constraints on the mass scales nor on the coupling constants. 



Chapter 5 
Conclusions 



In this thesis, gravitino production via different mechanisms has been discussed 
into details. We have first introduced the theories of inflation and of supergravity, 
which provide with the background for the research papers here included. These are 
studied in Chapters [3] and IH 

In the article PQ, the scattering of gauge bosons into gravitinos is investigated 
both in the gauge and mass eigenstates. If the W is massive, a term which grows 
with energy is produced in the cross section. This can accordingly violate the uni- 
tarity of the theory at large centre of mass scales. We show that this happens 
because the longitudinal degrees of freedom of the off-shell W, which is exchanged 
in the annihilation channel, are identically cancelled by the supergravity vertex. As 
proven in Section I3.1.2[ this implies that the on-shell W bosons become strongly 
interacting at the supersymmetry breaking scale. At this point, one could wonder 
whether missing diagrams might cancel the unexpected terms. We show however 
that inclusion of processes with the neutral Higgs doublet of the MSSM does not 
change the result. 

There are basically two perspectives from the above. Supergravity is an ef- 
fective theory, which is the limit of some unknown unified theory and is valid only 
under the SUSY breaking scale. Accordingly, our result can be used in phenomenol- 
ogy. Corrections to previous results (see (SH]) can be thus calculated. For example, 
the Boltzmann collision factor of the gravitinos which are produced in the thermal 
bath might be significantly modified. 

However, one might expect that the divergences would be cancelled by some 
yet unknown mechanism. For instance, additional diagrams might come from inter- 
action terms in the Lagrangian (12.661) which are usually not included. Moreover, the 
contribution of the hidden sector should be taken into account as well, since here 
both the EW gauge symmetry and supersymmetry are broken. This is a clear differ- 
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ence with respect to the Standard Model, where similar divergences in the scattering 
W^W~ — 7- W^W~ are cured uniquely by diagrams with the Higgs. 

In the paper [2] we instead discuss gravitino production from the decay of a 
scalar field 5*. At the end of inflation, 5* rolls down its potential and, when the mass 
ms becomes comparable to the Hubble rate, it starts coherent oscillations about 
the minimum of the potential. This mechanism generates decay into gravitinos. We 
first consider the dynamics of the scalar field. Since the Universe is expanding, the 
minimum of the potential is time-dependent. Accordingly, S tracks the minimum 
and a displacement 58 from the minimum is generated. This produces a very differ- 
ent oscillation amplitude with respect to previous analyses. We calculate the related 
corrections in a general setting, with a generic quadratic potential V{S). 

S is then embedded in a model of gauge mediation with metastable vacuum, 
that is provided with an explicit hidden sector of the O'Raifeartaigh type. The scalar 
field breaks supersymmetry via a linear term in the superpotential, and the Kahler 
potential for S", K{S), has a negative non-minimal term. We calculate the radiative 
and supergravity corrections to the scalar potential V{S), which come from one-loop 
diagrams of the scalars and fermions of the hidden and messenger sector. We show 
that they i) lift the fiat direction at the tree level by generating a mass term for S; 
ii) provide with an explicit form of the arbitrary cutoff A of the theory in terms of 
the parameters of the hidden sector; iii) induce a linear term in K{S) which makes 
the field deviate from the origin. 

By using the potential V{S) thus obtained, we calculate the yield variable 
of S and of the gravitino, by taking into account the correction factors previously 
obtained. By demanding that all the gravitinos thus produced constitute the main 
contribution to the amount of the observed Cold Dark Matter density, we find 
constraints on the reheating temperature Tr. A striking feature of the gravitino 
density and of the bounds fl4.102p and fl4.106p is their extreme sensitivity to the 
parameters of the hidden sector, especially to the coupling A of the O'Raifeartaigh 
model. This means that it is easy to account for gravitino dark matter with an 
arbitrarily low reheating temperature. 

These results suggest some directions of future research. In particular, the 
model might be completed with a source of baryogenesis, for example the Affleck- 
Dine mechanism. Also the impact of other sources of gravitino dark matter, for 
instance the decay of MSSM fiat directions, should be taken into account. Further 
constraints would then be imposed on a model that is a well-defined perturbative 
framework, by virtue of the explicit O'Raifeartaigh hidden sector. 

All the above perspectives are object of actual investigations. 
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Appendix A 

Notations and Feynman rules 



Here we provide with our notations and conventions, and with the Feynman rules 
of the MSSM and of supergravity which are used in Chapter [31 

A.l Notations and conventions 

Regarding our conventions, we follow basically the reviews by Haber and Kane [12] 
and Martin [H]. Throughout this thesis, the main subject is particle physics, even 
though we deal with cosmology. Therefore, it is a good approximation to adopt the 
flat space-time Lorentz metric 



with the standard convention of particle physics for the signature. This implies 
= m? for a particle with four-momentum p and mass m. The spinors in both 
the Standard Model and in supersymmetric Lagrangians are usually written first in 
two-component notation, since each Dirac fermion has left- and right-handed parts 
with completely different electroweak interactions. This way the degrees of freedom 
are treated differently already from the beginning. Accordingly, we introduce the 
Pauli matrices 



A two-component (left-handed) complex Weyl spinor transforms under actions 
M G SL(2, C), that is under the matrix representation of the Lorentz group S0(1, 3). 



= diag{l, -1,-1,-1) 



(A.l) 




(A.2) 
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The right-handed spinor transforms under M*, and similarly for and one 
has and (M~^)*. The spinor indices are raised and lowered according to the 
antisymmetric tensor e"^ = — e^": 

(A.3) 

where i]a is another Weyl spinor, since 

= _,21 ^ ^ ^ ^ . ^ ^ ^11 = ,22 ^ Q _ 

The fact that the Weyl spinors are anticommuting and the antisymmetric nature of 
e"'^, compensate each other in the following remarkable result: 

= Cea^V^ = -V^e^pC = V^e^aC = ■ (A.5) 
The Dirac equation in two-component notation can be written as 

{^.pT^ip = mff , {o^v^^^fif^ = , (A.6) 

where 

a^ = (l,a), ^^^ = (1,-^). (A.7) 
We can now rewrite Eq. (lA.6P in four-component notation, namely 

(7% -m)ij = 0, (A.8) 

by defining the spinor 

ta) ■ (A-9) 



The Dirac (or gamma) matrices in Eq. llA.SP are written in the Weyl basis as 

"o")- 7. = nV7V=(o' (A. 10) 

The above is called the chiral representation of the gamma matrices, and it is very 
useful in this context, where we want to separate the left- and right- polarizations 
of the Dirac spinor. To this aim, let us define the left- and right- handed projection 
operators: 

Pl,r = I {iTl') , (A.ll) 



so that is it possible to write 



where we have defined ipL,R = Pl,b.4'- 
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A.2 Feynman rules for the MSSM and SUGRA 



The relevant Feynman rules which are used in the scattering calculations of paper 
[1] are here summarised. The propagators in the R^ gauge and the vertices of the 
Minimal Supersymmetric Standard Model and of Supergravity follow Refs. 
and m\. 



In the following, Mp = (SttCtv)"^/^ = 2.43 x 10^^ GeV is the reduced Planck 
mass, g and /"^'^ are, respectively, the coupling constant and the structure constants 
of a generic gauge group. We point out that our expressions differ from those of 
Ref. [42] by an imaginary factor i, since the squared amplitude must be real. This 
is justified by the different normalization of the Lagrangians of Supergravity and of 
the MSSM. All the momenta are considered to run into the vertex and the fermion 
flow is denoted by a thin line with an arrow. 



A.2.1 Propagators and interaction vertices 



—I 



WAAAAAAr = ^2_^2 



6, u 



—I 



kp'—m'l 



v,u + - 1; 



knku 



— ^m? 



AA 



6^ 



ab 



3 I 

Figure A.l: The propagators of a vector boson with mass niA and of a scalar with 
mass m$. 



AAAAAA/VW = ^- , 



i^Ykp + mp) ^ij 



^2 _ 



Figure A.2: The propagators of a Majorana fermion with mass tum and of a Dirac 
fermion with mass m/). 
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Figure A. 3: Coupling of 3 gauge bosons. 





Figure A. 4: Gauge boson-gaugino-gaugino. 



H 



gmw cos - a)r]a^ ] 




wz 



Figure A. 5: Higgs-gauge bosons. 
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Figure A. 6: Gauge boson-chargino-neutralino. 
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Figure A. 7: Gauge boson-gaugino-gravitino. 
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Figure A. 8: Gauge boson-gauge boson-gaugino-gravitino. 
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Figure A. 9: Scalar- Weyl fermion-gravitino. 



